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REVISED  IMPACT  DYNAMIC  DESIGN-ANALYSIS  METHOD  (RIDDAM) 


INTRODUCTION 

BACKGROUND 

DP  AM 

A  particular  dynamic  design-analysis  method.  DDAM,  has  been 
applied  to  selected  equipment  on  Navy  ships  since  I960.  The  DDAM  is 
based  on  measurements  of  the  responses  of  shipboard  equipment  to 
controlled  tests  against  ships  with  underwater  explosions. 
Measurements  from  a  variety  of  ships  showed  that  the  responses  of 
equipment  to  the  shock  from  an  underwater-explosion  attack  could 
be  approximated  fairly  well  by  considering  onlv  the  location  of  the 
equipment  on  the  ship,  the  weight  of  the  equipment,  and  its  natural 
frequency  of  vibration.  Collected  measurements  were  fitted  to  simple 
formulas  depending  on  location,  weight,  and  frequency  that  could  be 
used  to  check  the  ability  of  equipment  to  resist  shocks  from 
underwater  explosions 

The  DDAM  can  be  looked  at  in  modern  terms  as  an  approximate 
method  of  conducting  a  substructure  analysis  of  a  complex  dynamic 
system  In  the  DDAM  all  of  the  effort  in  the  analvsis  is  directed 
.oward  one  particular  substructure  representing  the  equipment  of 
interest  The  main  structure  consisting  of  all  the  rest  of  the  ship,  is 
taken  as  a  standard  structure  whose  properties  do  not  need  to  be 
specially  determined  The  experimental  measurements  of  the  shock 
that  the  ship  can  transmit  to  equipment  of  different  weights  and 
frequencies  during  an  underwater-explosion  attack  take  the  place  of 
an  analysis  of  the  ship  s  structure 

The  effort  of  the  analvst  in  appiving  the  DDAM  is  limited  to 
determining  the  frequencies  and  weights  associated  with  one 
particular  item  of  equipment  The  formulas  from  the  DDAM  then 
provide  the  analyst  immediaieiv  with  the  responses  that  would  be 
expected  at  each  weight  and  frequency  if  the  equipment  were 
mounted  at  a  particular  location  on  a  standard  ship  during  an  attack 
against  the  ship  with  an  underwater  explosion  The  analyst  performs. 
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in  effect,  a  substructure  analysis  without  having  to  analyze  the  main 
structure  or  consider  location  of  the  equipment  beyond  the  broad 
categories  given  in  the  formulas 

1DDAM 

Calculations  for  the  DDAM  are  based  on  the  assumption  that  the 
shock  from  an  underwater  explosion  produces  a  translation  of  the  ship 
in  a  single  direction,  with  all  of  the  supports  for  the  equipment 
moving  identically.  This  is  generally  a  good  approximation  for 
compact  equipment  or  for  explosions  that  take  place  in  the  water 
some  distance  from  the  ship 

An  extension  of  the  DDAM.  referred  to  as  IDDAM  (Impact 
DDAM)  was  developed  in  1977  in  response  to  a  request  for  a  dynamic 
design-analysis  method  that  could  be  applied  to  the  case  of  the  impact 
of  exercise  (nonexploding)  torpedoes  against  the  hull  of  a  submarine 
Here  the  shock  applied  to  the  several  supports  for  a  particular  item  of 
equipment  attached  to  the  hull  could  vary  greatly,  from  a  severe 
shock  to  a  support  directly  in  way  of  the  impact  to  shocks  of  smaller 
and  smaller  severity  with  increase  in  distance  from  the  point  of 
impact. 

Determining  the  response  of  a  structure  to  independent  motions 
of  multiple  supports  involved  only  a  minor  elaboration  of  DDAM 
However,  the  weight  of  the  structure  had  to  be  distributed  among  the 
separate  supports  in  a  self -consistent  fashion  in  order  to  allow  the 
simple  formulas  involving  weight  and  frequency  to  be  applied.  When 
all  of  the  supports  move  identically,  as  in  DDAM.  the  full  weight  reacts 
against  the  common  motion,  but  when  the  motions  are  different  the 
weight  appears  as  a  combination  of  direct  and  cross  terms  among  the 
supports,  adding  appreciably  to  the  complexity  of  the  method. 

REVISION  OF  THE  IDDAM 
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RIDDAM 

The  present  report  describes  a  new  version  of  the  IDDAM 
referred  to  here  as  RIDDAM  'Revised  IDDAM I.  The  RIDDAM  simplifies 
the  method  of  calculating  the  effective  weight  of  a  structure  with 
multiple  supports. 
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PROCEDURE 


Mathematical  model 

The  equipment  (substructure)  is  modeled  as  a  linear,  elastic, 
and  undamped  structure.  Details  of  the  model  depend  on 
requirements  of  the  computer  program  that  is  used. 

Frequencies 

Frequencies  are  found  by  doing  a  standard  normal-mode 
analysis  of  the  model  with  fixed  interfaces  at  the  points  of  attachment 
to  the  ship. 

Participation  factors 

The  normal-mode  analysis  provides  a  set  of  factors  showing 
how  the  motion  of  each  support  participates  in  producing  a  response 
of  each  of  the  normal  modes. 

Modal  weights 

The  critical  feature  of  both  the  DDAM  and  the  IDDAM  is  the 
assignment  of  an  effective  mass  to  each  normal  mode  of  vibration 
The  effective  mass  is  used  as  a  measure  of  the  reactive  force  that  the 
mode  applies  to  the  supports  for  the  equipment  when  the  equipment 
is  accelerating  in  response  to  a  shock.  In  the  RIDDAM  this  effective 
mass  is  obtained  simply  by  adding  mass-normalized  participation 
factors  by  absolute  value  and  squaring  the  sum. 

Irmy.tf 


As  for  the  DDAM,  the  inputs  for  the  IDDAM  are  formulas 
showing  the  response  of  a  normal  mode  of  vibration  as  a  function  of 
its  weight,  frequency,  and  location  on  board  ship  In  the  DDAM,  all  of 
the  supports  are  assumed  to  move  together  in  one  particular  direction. 
The  IDDAM  relaxes  this  assumption  by  providing  different  inputs  for 
different  supports  The  RIDDAM  uses  the  same  inputs  as  the  IDDAM 


Computer-aided  analysis 

The  original  DDAM  was  developed  at  a  time  when  structural- 
analysis  computations  were  regularly  done  using  pencil,  paper,  and  a 
desk  calculator.  Its  procedures  have  been  carried  over  to  the  present 
IDDAM. 

First,  turn  on  your  computer.  In  the  following,  the  equations  for 
the  IDDAM  have  been  written  using  the  matrix  notation  that  is 
common  for  computer-aided  analysis,  on  the  assumption  that  a 
computer  rather  than  a  desk  calculator  will  be  used  to  analyze  the 
structure.  Appendix  A  gives  some  of  the  elements  of  matrix  algebra 
for  readers  who  are  not  entirely  familiar  with  matrix  notation. 
Appendix  B  lists  a  computer  program  that  can  be  used  to  analyze 
simple  structures  on  any  computer  that  has  a  compiler  or  interpreter 
for  the  BASIC  computer  language. 

The  rewritten  equations  for  the  IDDAM  produce  a  more  compact 
notation  and  simpler  calculations  but  make  no  change  in  its 
procedures. 

Effective  weight 

The  RIDDAM  does  make  a  substantial  change  in  procedures  bv 
the  method  it  uses  to  determine  effective  weight,  however 

In  both  the  DDAM  and  the  IDDAM  the  weight  of  a  normal  mode 
of  vibration  is  considered  as  reacting  against  supports  only  to  the 
extent  that  the  supports  are  driven  by  a  shock  input  The  weight  of  a 
mode  thus  varies  with  direction  of  the  input  in  the  DDAM  and  varies 
with  the  distribution  of  shock  severities  among  the  supports  in  the 
IDDAM. 

The  method  of  determining  the  weight  of  a  mode  in  the 
RIDDAM  presumes  that  the  mode  reacts  against  all  of  the  supports 
whenever  it  responds  to  a  shock  delivered  through  any  of  them  The 
effective  modal  weight  in  the  RIDDAM  is  thus  independent  of  the 
particular  shock  that  is  applied.  Calculations  are  not  only  simplified 
but  reduced  to  the  more-rational  basis  of  assuming  that  the  structure 
reacts  against  all  of  its  supports  whenever  it  is  forced  to  undergo 
accelerations  from  anv  cause 
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Responses  of  the  modes  are  combined  to  estimate  forces, 
stresses,  strains,  and  deformations  throughout  the  model  of  the 
equipment.  The  stresses  and  deformations  at  critical  locations  are 
compared  with  allowable  values  to  assess  the  ability  of  the  equipment 
to  withstand  a  particular  shock.  Deformations  produced  by 
differences  in  the  motions  of  redundant  supports  are  included  in  both 
the  IDDAM  and  the  RIDDAM. 

OUTLINE 

The  report  begins  with  a  brief  overview  of  the  theoretical  basis 
for  the  normal-mode  method  used  in  the  RIDDAM.  It  continues  with  a 
discussion  of  some  of  the  practical  compromises  that  are  necessary  to 
keep  its  implementation  simple,  and  a  description  of  some  of  the 
problems  that  the  compromises  generate.  There  are  some  worked-out 
examples  of  analyses  of  simple  structures  by  the  RIDDAM.  The  object 
is  to  supply  enough  information  so  that  an  analyst  can  apply  the 
RIDDAM  to  his  own  item  of  equipment  by  following  the  discussions 
and  examples  given  here 
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THEORY 


MATHEMATICAL  MODEL 

Balance  of  forces 

Most  computer  programs  for  dynamic  analvsis  represent  the 
structure  by  a  matrix  equation  of  the  form 

Ma  +  K  x  -  Fu.  Hi 

where  ait)  is  a  time-varying  column  vector  of  accelerations  at  selected 
points  (nodes)  on  the  structure  and  lit)  is  a  column  vecto.r  of 
displacements.  The  symmetric  matrices  M  and  K  represent  the  mass 
and  stiffness  of  the  structure.  The  column  vector  u(t)  consists  of 
displacements  of  the  supports  for  the  structure;  the  rectangular  matrix 
of  support  stiffnesses  F  converts  displacements  of  the  supports  into 
forces  applied  to  the  coordinates  of  the  structure  The  equation 
balances  these  forces  against  the  forces  Ma  and  Kx  from  the  mass  and 
stiffness  of  the  structure. 

The  motions  represented  by  Equation  1  can  be  translations  in 
any  direction  or  rotations  about  any  axis,  with  elements  of  M.  K.  and  F 
being  mass,  mass  moment  of  inertia,  force  per  deflection,  or  moment 
per  angle.  It  is  assumed  that  all  the  translations  and  rotations  are 
small,  the  stiffnesses  are  constant,  and  there  is  no  damping  in  the 
structure. 

Discussions  here  are  framed  in  terms  of  translations 
(displacements  and  forces)  with  the  understanding  that  rotational 
motions  (angles  and  moments)  are  also  included 

OtheL-Lecresemations 

Some  methods  of  structural  analysis  are  based  on  other  than 
the  force-balance  equation.  For  example,  the  displacement-balance 
equation 

K-l  M  a  ♦  x  =  K"1  F  u  <2 ' 
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be  a  diagonal  matrix.  There  are  many  different  numerical  methods 
for  generating  a  Q  that  will  satisfy  Equations  3  and  4  simultaneously 
Which  is  best  depends  on  the  size  of  the  problem,  the  distribution  of 
the  nonzero  elements  in  M  and  K,  and  how  many  elements  of  Q  are  to 
be  calculated.  Many  programs  offer  a  choice  of  methods  and  may  also 
have  ways  of  decreasing  the  sizes  of  the  structural  matrices  to  speed 
up  the  calculation. 


If  the  matrix  0  is  used  to  transform  the  accelerations  and 
displacements  by 


a  -  Q  b  i  -  0  y. 

and  if  Equation  1  is  premuitiphed  by  Q^.  there  results 

b  ♦  p2  y  _  qT  p  u  1 6' 

Because  P2  is  diagonal,  each  element  of  the  acceleration  blu  and 
displacement  y(t)  on  the  left  side  represents  an  independent  oscillator 
with  a  fixed-base  angular  frequency  given  by  the  square  root  of  its 
element  of  P2  The  responses  of  the  oscillators  can  be  converted  back 
to  motions  of  the  structure  bv  Equation  5 


Equation  6  can  be  reduced  to  standard  form  for  each  oscillator 
by  writing  it  as 


b  +  P2  y  =  P2  RT  u,  i ,7 1 

where 

rT  .  P'2  QT  p  .  8  i 

is  a  rectangular  matrix  called  the  participation  factor.  The  matrix  RT 
converts  the  displacements  ui.U  of  the  supports  to  a  generalized 
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has  been  popular  for  use  with  a  desk  calculator.  Here  (the  inverse 
of  the  stiffness  matrix)  is  a  matrix  of  influence  coefficients  that  shows 
the  displacements  of  all  the  points  on  the  structure  when  unit  force  is 
applied  to  one  of  the  points.  The  support  influence  matrix  K_1F 
represents  the  position  of  the  structure  produced  by  the  motion  of  its 
supports.  The  equation  balances  these  equilibrium  displacements 
against  the  dynamic  displacements  that  include  forces  from  the 
accelerations  of  the  masses. 

The  displacement-balance  equation  is  not  well  suited  to 
computer-aided  analyses  because  the  matrix  K_IM  is  full  of  nonzero 
elements  and  is  not  symmetric.  A  computer  needs  less  memory  and 
can  run  faster  using  the  matrices  M  and  K.  both  of  which  are  sparse 
and  symmetric. 

Discussions  here  are  in  terms  of  the  force-balance  equation  as 
the  equation  most  common  in  computer  programs  for  dynamic 
structural  analysis. 

NORMAL  MODES  OF  VIBRATION 


Many  computer  programs  begin  a  normal-mode  analysis  by- 
generating  a  square  or  rectangular  matrix  Q  that  makes 

0T  M  0  -  I  1 3  • 

be  a  unit  matrix.  Normalizing  on  the  mass  matrix  in  tins  wav- 
produces  a  simple  and  efficient  analysis,  and  has  been  assumed  in  the 
following  discussion.  A  later  paragraph  will  consider  the  alterations 
that  are  necessary  when  other  methods  of  normalization  are  used 


Diagonalizing  the  stiffness  matrix 


The  matrix  Q  can  also  be  chosen  to  make 


displacement  representing  the  motion  of  the  base  of  each  one  of  the 
normal-mode  oscillators.  The  participation-factor  matrix  has  as  many 
columns  as  there  are  separate  supports  for  the  structure,  and  as  many 
rows  as  there  are  normal  modes. 

Each  of  the  normal  modes  in  Equation  7  responds  as  if  it  were  a 
single  oscillator  supported  from  a  base  whose  displacement  is  given 
by  one  of  the  elements  of  RTu(t).  That  is,  if  y(t)  is  the  displacement 
and  b(t)  is  the  acceleration  of  one  of  the  normal  modes,  the  equation 
for  its  response  is 

b  ♦  p2  v  =  p2  j-T  u.  1 9 1 

where  rT  is  the  row  from  RT  that  corresponds  to  the  elements  y«t)  and 
bit),  and  P2  is  the  square  of  the  angular  frequency  of  the  mode 

MODAL  MASS 

What  the  modal  mass  does 

A  force  is  required  to  produce  an  acceleration  of  a  real  structure 
that  has  mass.  This  force  reacts  against  the  supports  for  the  structure 
and,  in  all  practical  cases  where  the  supports  are  attached  to  some 
other  structure,  produces  motions  of  the  supporting  structure 

The  displacements  ult)  of  the  supports  need  to  be  considered  as 
consisting  of  two  components  One  component  is  the  motion  'or  lack  of 
motion)  that  would  be  expected  if  the  structure  were  not  in  place  The 
other  component  is  the  motion  produced  by  the  reaction  force  as  the 
structure  accelerates. 

If  the  structure  is  analyzed  into  its  normal  modes,  each  normal 
mode  needs  to  be  assigned  a  mass  so  its  reaction  force  can  be 
determined. 


A  simple  example  of  modal  mass 

Suppose  a  single  mass  VI  is  supported  by  a  set  of  parallel 
springs  having  a  total  stiffness  k  The  equation  of  motion  can  be 
written 


where  fT  is  a  row  vector  whose  elements  are  the  stiff  nesses  of  the 
individual  springs  and  u(t)  is  a  column  vector  of  the  displacements  at 
the  points  where  each  spring  is  attached  to  some  other  structure 

Follow  the  method  of  analysis  given  in  the  preceding  sections  as 
if  the  mass  M  and  the  stiffness  K  were  matrices  with  just  one  row  and 
one  column.  First  define  a  modeshape  Q  that  makes 

QMQ  ■  1.  '11' 

Notice  that  Q  is  the  inverse  square  root  of  M.  and  that 

QKQ  -  K/M  *  P2  1  12' 

is  the  square  of  the  natural  frequencv  of  the  mass  as  supported  by  all 
the  springs. 

Continue  to  find  the  participation  factor 

rT  -  p-2  Q  fT  i  1  3  i 

and  reduce  it  to 

rT  =  (K/M)’1  M'1/2  fT  -  vM  fT/k.  <M' 

Summing  the  elements  of  the  participation -factor  vector  gives  the 
square  root  of  M.  and  squaring  the  sum  recovers  M  as  the  mass  of  the 
single  mode  of  the  structure 

The  same  procedure  can  be  applied  to  Equation  9  to  determine 
an  effective  mass  for  the  mode  whose  acceleration  is  b(t)  That  is  add 
the  elements  of  the  row  from  the  participation-factor  matrix  and 
square  the  sum  to  find  the  mass  that  must  be  accelerated  when  the 
structure  responds  with  motions  in  one  particular  normal  mode  of 
vibration. 
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that  with  independent  motions  of  multiple  supports,  the  reactions 
occur  and  the  supports  are  moved  whenever  a  mode  is  accelerated, 
whether  the  acceleration  was  initiated  by  a  shock  delivered  through  a 
particular  support  or  not. 

Using  the  square  of  the  sum  of  the  absolute  values  of  the  mass- 
normalized  participation  factors  as  a  measure  of  the  mass  of  a  normal 
mode  of  vibration  is  the  only  substantial  change  that  the  RIDDAM 
makes  in  the  IDDAM.  The  IDDAM  distributed  the  mass  of  each  mode 
among  the  supports  in  accordance  with  the  shock  severity  at  each 
support.  In  the  RIDDAM  each  mode  has  a  mass  that  is  independent  of 
the  shocks  that  are  applied 


In  the  original  DDAM.  the  shock  motion  of  a  ship  was  considered 
to  be  a  rigid-body  translation  in  a  single  direction.  All  the  supports 
for  an  item  of  equipment  underwent  the  same  motion  and  the  masses 
of  the  modes  were  the  effective  masses  reacting  against  that  common 
motion.  In  the  RIDDAM  each  support  is  taken  as  attached  to  a 
structure  that  can  move  independently  of  the  other  supports  and  the 
effective  mass  of  each  mode  reacts  against  ail  of  the  supports  in 
whatever  direction  applies  at  each  support. 

It  is  recognized,  of  course,  that  points  located  near  one  another 
on  the  structure  of  a  ship  are  neither  tied  rigidly  together  IDDAM)  nor 
completely  independent  (RIDDAM)  The  object  of  the  procedures 
adopted  in  each  case  is  to  provide  a  reasonable  and  simple  method  of 
estimating  the  effects  of  the  reaction  forces  without  having  to  extend 
the  analysis  into  the  details  of  the  supporting  structure. 


The  example  shown  in  Equations  10  to  14  was  for  a  single  mass 
moving  in  a  single  direction  Later  examples  will  show  that  the 
procedure  of  adding  influence  coefficients  by  absolute  value  and 
squaring  the  sum  gives  an  effective  mass  for  a  mode  even  when  a 
structural  model  includes  translations  in  different  directions  or 
rotations  about  different  axes. 
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The  advance  of  the  DDAM  was  in  pointing  out  that  the  effect  of 
reaction  forces  on  a  particular  structure  can  be  found  experimentally 
by  mounting  substructures  to  it  and  measuring  the  responses  of 
substructures  having  different  masses  and  frequencies.  The  results 
from  the  tests  can  then  be  applied  to  predict  the  responses  of  modes 
of  other  items  mounted  to  the  structure  provided  their  masses  and 
frequencies  are  known.  The  supporting  structure  thus  does  not  need 
to  be  analyzed  in  order  to  estimate  the  shock  response  of  items 
attached  to  it. 

Some  approximation  is  inherent  in  this  assumption,  especially  if 
responses  measured  for  one  particular  structure  are  applied  to 
different  but  similar  structures.  Supporting  structures  can  differ 
because  of  other  substructures  mounted  near  the  substructure  of 
interest,  or  even  as  a  result  of  other  modes  of  the  substructure  being 
analyzed.  However,  the  advantage  of  being  able  to  estimate  the 
response  of  a  particular  item  of  equipment  without  having  to  make  a 
detailed  analysis  of  the  supporting  ship  s  structure  makes  the  method 
attractive  despite  the  potential  errors. 


Squaring  the  sum  of  the  elements  in  a  row  of  the  influence- 
coefficient  matrix,  as  shown  in  the  simple  example  above,  will  give  an 
effective  mass  for  the  mode  if  all  the  reactions  are  in  the  same 
direction,  as  they  were  for  the  parallel  springs  in  the  example  If 
some  of  the  reactions  are  positive  and  some  are  negative,  the  squared 
sum  will  give  an  effective  mass  producing  a  net  force  on  all  of  the 
supports  taken  together  This  is  the  effective  mass  used  m  the  DDAM. 
where  differences  in  the  motions  of  different  supports  are  not 
considered. 

If  each  point  of  support  is  taken  as  attached  to  a  separate 
structure,  the  reactions  will  produce  independent  motions  of  each 
support  in  accordance  with  the  direction  and  magnitude  of  each  force 
In  this  case  the  effective  mass  acting  against  all  the  supports  should 
be  determined  by  adding  the  participation  factors  by  absolute  value 
(that  is,  ignoring  their  directions)  and  then  squaring  the  sum  Notice 
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The  IDDAM  inputs  represent  linear  accelerations  as  a  function 
of  effective  mass  and  frequency.  Applied  through  different  supports, 
these  accelerations  can  induce  rotational  and  cross-axis  responses  of  a 
structure.  However,  no  inputs  are  supplied  to  show  how  the  response 
of  a  structure  to  a  rotation  of  one  of  its  supports  might  depend  on  the 
mass  moment  of  inertia  about  that  support. 

The  RIDDAM  uses  the  inputs  developed  for  the  IDDAM.  As  a 
result  of  the  lack  of  rotational  inputs,  it  is  necessary  that  the  structure 
be  modeled  so  that  its  response  depends  on  translational  inputs  only. 
This  can  be  done  in  either  of  two  ways:  A  support  can  be  taken  as 
irrotational,  so  that  moments  applied  to  it  do  not  cause  any  rotation. 
Alternately,  the  joint  at  the  support  can  be  taken  as  pinned,  so  that 
the  structure  can  rotate  freely  about  that  point  without  producing  any 
rotation  of  the  supporting  structure. 

INPUTS 

Differential  equations 

The  inputs  to  the  analysis  consist  of  predetermined  solutions  to 
the  equations 


Me  ♦  M  P2  z  -  M  P2  u,  1 1  Vi 

where  dt)  is  the  acceleration  and  z(t)  is  the  displacement  of  an 
oscillator  having  mass  M  and  frequence'  P.  Here  u< t i  is  the 
displacement  of  the  support  for  the  oscillator,  including  both  the 
motion  induced  by  an  external  force  and  the  motion  resulting  from  the 
reaction  of  the  mass  of  the  oscillator  against  its  support. 

Because  of  the  reaction,  the  acceleration  of  the  oscillator  varies 
with  its  mass  as  well  as  with  its  frequency  In  the  DDAM.  all  the 
support  motions  are  taken  as  identical,  so  only  one  function  is  needed 
for  the  base  motion  ult),  but  in  the  IDDAM  the  supports  can  move 
independently  and  there  is  a  separate  Equation  15  for  the  effect  of 
each  support  on  each  mode. 
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The  equation  for  the  response  of  one  of  the  normal  modes, 

b  +  p2  v  -  p2  rT  u,  ( 


uses  a  row  from  the  influence-coefficient  matrix  to  take  a  linear 
combination  over  the  motions  of  all  of  the  supports  and  uses  that 
combination  to  find  the  response  of  the  mode.  The  process  can  be 
reversed  by  calculating  responses  first  and  then  taking  the  linear 
combination  to  obtain 


where  bit)  is  the  acceleration  of  the  mode  and  c(M,  P,  t)  is  a  column 
vector  of  the  solutions  to  Equation  15  for  the  mass  and  frequencv  of 
the  mode  and  for  the  shock  applied  to  each  of  the  supports. 

If  the  response  accelerations  c(M,  P,  t)  for  each  support  are 
already  provided.  Equation  17  will  give  the  responses  for  the  modes  of 
the  substructure  without  the  necessity  of  solving  any  differential 
equations. 

ACCELERATIONS.  FORCES.  STRESSES.  AND  DEFORMATIONS 


Scale  the  matrix  of  modeshapes  to 


where  Bit)  is  a  diagonal  matrix  of  the  accelerations  bit)  of  the 
individual  modes  and  Ail)  is  now,  by  Equation  5.  a  matrix  whose 
columns  represent  the  accelerations  of  the  coordinates  of  the  structure 
in  each  mode.  It  is  convenient  at  this  point  to  keep  the  responses  in 
the  different  modes  separate  because  of  approximations  that  will  be 
made  later  in  combining  responses  in  different  modes. 


Forces 


Multiply  the  matrix  of  accelerations  by  the  mass  matrix  M  to 


obtain 

G  -  M  A  (19) 

as  a  matrix  whose  columns  represent  the  forces  applied  to  the 
structural  nodes  by  the  accelerations  of  each  mode 


The  forces  will  produce  time-varying  stresses  and  deformations 
in  the  structure  that  are  proportional  to  the  combination  of  forces  in 
each  mode.  The  stresses  and  deformations  at  critical  points  on  the 
structure  can  be  used,  to  estimate  the  ability  of  the  structure  to 
survive  the  accelerations  given  by  the  inputs. 

WARPAGE 

Additional  deformations 

The  displacement-balance  equation  (Equation  2). 

x  -  K-l  F  u  -  K-l  M  a  (20) 

shows  that  the  response  i(t)  of  the  structure  depends  on  the  motions 
of  its  supports  u(t)  as  well  as  on  its  acceleration  ait)  The  preceding 
modal  analysis  accounts  for  the  forces,  stresses,  and  deformations 
produced  by  acceleration  only.  It  is  necessary  to  supplement  the 
acceleration-generated  responses  with  a  check  of  any  additional 
deformations  <and  resulting  stresses)  that  may  be  produced  bv 
differences  in  motions  of  individual  supports 

Note  that  if  the  motions  of  the  supports  are  all  identical  'as 
assumed  in  the  DDAM)  there  will  a  displacement  response  but  no 
additional  deformation  of  the  structure.  There  also  is  no  warping 
deformation  if  the  supports  are  not  redundant  or  if  the  motions  of  the 
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supports  can  be  described  bv  a  combination  of  rigid-body 
displacements  and  rotations 


The  structure  will  resist  warping  deformation  by  generating 
static  forces  that  tend  to  decrease  the  differences  between  the  motions 
of  its  supports.  Again,  as  in  the  case  of  the  effect  of  modal  mass  on 
acceleration  responses,  an  estimate  of  the  effect  of  the  reaction  forces 
can  be  made  in  terms  of  an  experimental  or  analytical  decrease  in 
warping  with  increasing  static  reactions  of  a  structure  against  relative 
motions  of  a  standard  or  typical  supporting  structure. 

SUMMARY 

The  preceding  theory  shows  how  the  responses  of  substructures 
can  be  estimated  under  standard  conditions  while  mounted  to  a 
particular  supporting  structure  The  characteristics  of  the  supporting 
structure  are  evaluated  by  determining  the  responses  it  can  produce 
in  oscillators  having  different  masses  and  frequencies  and  the  static 
warpage  it  can  produce  in  redundant  supports  with  different 
reactions. 

The  structure  of  interest  is  then  analyzed  mto  oscillators  having 
different  masses  and  frequencies  and  into  reactions  against  relative 
motions  of  its  supports  Direct  comparison  with  the  responses 
determined  for  the  supporting  structure  then  allows  responses  of  the 
structure  to  be  estimated 

The  theory  has  the  potential  for  being  exact  if  the  supporting 
structure  were  identical  from  one  substructure  to  the  next.  In 
application,  however,  supporting  structures  lend  to  be  similar  but  not 
identical.  Then  the  accuracy  of  the  method  depends  on  how  much 
variety  is  allowed  in  supporting  structures  before  a  new  set  of 
characteristics  must  be  applied 
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APPLICATION 


UNITS 

Convenience 

The  IDDAM  is  based  on  user-friendly  engineering  units  in  which 
Newton's  law  is  satisfied  for  forces  in  pounds,  masses  in  pounds,  and 
accelerations  in  multiples  of  the  acceleration  of  gravity.  The 
acceleration  of  gravity  is  defined  as  356  inches  per  second  per  second 
(9,8044  meters  per  second  per  second  with  one  inch  equal  to  0.0234 
meters).  A  pound  of  force  is  then  4.447  Newtons,  while  a  pound  of 
mass  is  0.4336  kilograms. 

Inputs 

The  tabulated  inputs  for  the  IDDAM  give  accelerations  m 
multiples  of  the  acceleration  of  gravitv  (g)  as  a  function  of  weights  in 
pounds  and  frequencies  in  Hertz  'cycles  per  second). 


Most  computer  programs  will  not  allow  forces  and  masses  both 
to  be  entered  in  pounds  A  common  convention  for  using  such 
programs  in  the  inch-pound-second  system  is  to  enter  masses  m  units 
of  356  pounds.  That  is.  weights  in  pounds  are  divided  by  356  to 
obtain  masses  for  entry  into  the  computer  program  Forces  can  then 
be  entered  in  pounds,  with  stiffnesses  in  pounds  per  inch  and  elastic 
moduli  in  pounds  per  square  inch  There  is  no  accepted  name  for  the 
unit  of  mass  that  weighs  356  pounds,  although  the  term  slinch  has 
been  suggested  as  a  parallel  to  the  term  slug'  that  is  used  for  the 
similar  unit  in  the  foot-pound-second  svstem 
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SCALING  OF  MODESHAPES 


Some  computer  programs  apply  scale  factors  to  the  modeshapes 
to  produce  a  modal  matrix  QS,  where  S  is  a  diagonal  matrix  of  scale 
factors.  One  popular  method  of  scaling  makes  the  largest  element  in 
each  column  of  QS  have  the  value  unity,  for  example.  The  theory 
shown  here  was  based  on  mass-normalized  modeshapes  in  which 
QtMQ  -  I  was  a  unit  matrix.  Modifications  are  needed  if  modeshapes 
have  been  scaled  otherwise. 


With  scaled  modeshapes. 

(QS )T  M  QS  -  S2  (211 

is  diagonal  but  not  a  unit  matrix.  Also, 

(QS)t  K  QS  =  S2  P2  122) 

is  diagonal  but  its  elements  are  no  longer  the  squares  of  the  natural 
frequencies.  Using  the  scaled  modeshapes  in  the  transformations 


a  -  QS  b 


i  -  QS  y 


converts  the  force-balance  equation 


M  a  •  Ki  =*  F  u 


S2  b  *  S2  P2y  -  'QS)T  F  u 


as  the  analog  of  the  normal-mode  equations. 


The  elements  of  S2  in  Equation  25  have  been  referred  to  as 
masses  of  the  normal  modes  and  even  used  as  if  they  were  masses  in 
some  calculations.  One  popular  computer  program  refers  to  them  as 
"generalized  masses'"  in  its  printed  output.  The  elements  of  S2  have 
completely  arbitrary  positive  values  that  depend  only  on  what 
criterion  was  used  to  scale  the  modeshapes  and  have  no  physical 
significance. 


The  scale  factors  contained  in  S  must  simply  be  removed  step- 
by-step  as  the  analysis  proceeds.  First,  divide  by  the  elements  of  S2 
to  obtain  the  squares  of  the  natural  frequencies  as 

p2  .  S'2  (S2  P2).  (26  ) 

Next,  recognize  two  kinds  of  participation  factors.  The  mass- 
normalized  participation  factor  must  be  recovered  from 

rT  .  s-1  IP'2  (OS  )T  F],  (27 ! 

The  responses  of  the  scaled  modes,  however,  must  be  found  from  the 
equation 

b  *  P2y  -  P2  [P-2  S'2  (QS)TFI  u,  <28 ) 

so  that  the  participation  factor  for  this  operation  must  be  obtained 
from 

r,T  »  P'2  S'2  (QS)T  F  129  i 

Notice  that  any  computer  program  that  displays  natural 
frequencies  must  have  calculated  a  value  of  S2  in  order  to  evaluate 
Equation  26,  so  it  should  not  be  necessarv  to  recalculate  the  scale 
factors  (Equation  2 1 )  in  order  to  find  what  factors  to  remove. 
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Avoiding  the  scaling 


Many  computer  programs  begin  a  normal-mode  analysis  by 
choosing  a  modeshape  matrix  that  reduces  QTMQ  to  a  unit  matrix 
With  such  programs  the  mass-normalized  modeshapes  are 
fundamental  and  subsequent  scaling  may  be  an  option  that  need  not 
be  selected. 

If  only  scaled  modeshapes  are  available  from  a  particular 
program,  it  may  be  more  efficient  to  remove  the  scale  factors  directly 
from  the  modeshapes  by  dividing  each  column  of  QS  by  its  element  of 
S  and  then  proceeding  m  a  mass-normalized  fashion,  rather  than 
dividing  repeatedly  by  the  scale  factors  as  results  are  being  calculated 

SUPPRESSION  OF  TIME  HISTORIES 


Complication 

The  theory  indicates  that  a  supporting  structure  must  be 
characterized  by  a  set  of  time -history  accelerations  c(M.  P.  t)  that 
differ  for  each  mass  M  and  frequency  P  of  oscillator  that  may  be 
supported  from  it.  Moreover,  the  time  histories  cannot  be  oscillations 
at  the  frequency  P,  because  P  is  a  fixed -base  frequency,  and  the 
essential  feature  of  the  supporting  structure  is  that  it  cannot  be 
considered  to  be  a  fixed  base. 

In  addition,  the  response  of  a  multimode  structure  with 
multiple  independent  supports  must  be  determined  by  superposing 
separate  nme-historv  responses  for  each  mode  and  each  support 
Such  combinations  were  impractical  for  routine  analyses  when  DDAM 
was  under  development  about  1960. 


The  IDDAM  follows  the  lead  established  by  the  DDAM  in 
suppressing  the  time  variation  of  the  response  acceleration  It 
characterizes  supporting  structures  in  terms  of  the  peak  'largest 
absolute  value]  of  acceleration  that  they  can  produce  lor  an  oscillator 
of  mass  M  and  fixed-base  frequency  P. 
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Suppressing  the  time  histories  allows  supporting  structures  to 
be  described  by  double-entry  tables  showing  peak  acceleration  as  a 
function  of  mass  and  frequency,  and  eliminates  the  superposition  of 
time  histories  for  multiple  modes  and  multiple  supports  in 
determining  responses.  It  leads  to  problems,  however,  in  estimating 
the  combined  effects  of  several  modes  whose  individual  responses  are 
given  only  in  terms  of  their  peak  values, 

TIME  HISTORIES  AND  COMBINATION  RULES 

Upper  and  lower  bounds 

It  is  clear  that  if  the  peak  responses  of  several  modes  were  all 
to  occur  in  the  same  direction  at  the  same  instant  of  time,  the 
combined  response  would  be  the  sum  of  the  individual  peaks  and  that 
this  sum  would  be  an  upper  bound  for  a  time-historv  combination  of 
responses.  A  lower  bound  may  also  exist  if  one  of  the  peak  values  is 
large  enough  to  dominate  the  combination  even  if  all  the  other  peaks 
happened  to  occur  simultaneously  in  the  opposite  direction.  Either 
occurrence  is  possible  but  unlikely  in  a  sum  of  time  histories. 

Combination  rules 

Simple  and  arbitrary  rules  can  be  defined  to  estimate  a  peak 
from  a  combination  of  responses  that  are  described  only  by  their 
individual  peak  values.  Any  rule  adopted  ought  to  define  a  combined 
peak  falling  in  the  range  from  the  lower  to  the  upper  bound.  The  rule 
should  select  a  likely  value  in  that  range  and  should  be  adjusted  to 
avoid  too  much  overconservatism  or  underconservatism. 

The  exact  form  of  a  combination  rule  is  not  important  in  areas 
of  a  structure  where  one  mode  predominates,  so  that  the  upper  and 
lower  bounds  do  not  differ  greatly  A  combination  of  engineering 
judgment,  common  sense  and  experimental  data  are  needed  to 
determine  combination  rules  when  several  modes  have  comparable 
peak  values  so  that  there  is  a  large  spread  between  upper  and  lower 
bounds.  Under  this  condition  the  best  rule  would  be  one  that  avoids 
obvious  conditions  of  great  overestimates  or  underestimates,  rather 
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than  a  rule  that  attempts  to  recreate  an  actual  combination  of  time- 
history  responses. 

COMBINING  PEAK  RESPONSES  FROM  DIFFERENT  SUPPORTS 


A  suitable  rule  for  estimating  the  peak  acceleration  b  of  one 
mode  from  a  collection  of  peak  response  accelerations  c  at  different 
supports  is 


where  r?  is  the  row  for  the  mode  from  the  participation-factor  matrix. 


Discussion 

Equation  30  is  identical  to  Equation  17  that  was  used  for 
combining  time  histories.  It  presumes  that  the  peak  responses  of  the 
mode  produced  by  the  motions  of  each  support  all  occur  at  the  same 
instant  of  time  and  in  directions  given  by  the  signs  of  the  elements  of 
the  participation  factor.  Such  a  combination  could  actually  occur 
under  torpedo  impact  if  the  torpedo  w'ere  to  strike  at  the  center  of  a 
symmetric  array  of  supports,  producing  large  responses  to  symmetric 
modes  of  the  structure  and  no  response  for  antisvm metric  modes 

Equation  30  also  reconciles  the  RIDDAM  with  the  DDAM  by 
matching  us  results  with  those  of  the  DDAM  for  conditions  in  which  all 
of  the  supports  would  have  identical  motions  in  a  single  direction. 

In  other  cases  the  rule  provides  an  estimated  peak  that  falls 
between  the  upper  and  lower  bounds  for  the  combination. 

COMBINING  PEAK  RESPONSES  FOR  MODES  OF  DIFFERENT 
FREQUENCIES 


An  estimate  of  the  peak  stress  or  deformation  from  a 
combination  of  peak  stresses  or  deformations  from  modes  of  different 


lrequencies  is  obtained  by  selecting  the  response  having  the  largest 
peak  in  absolute  value  and  adding  to  it  the  square  root  of  the  sum  of 
the  squares  of  the  peak  values  of  the  other  responses.  This  method  of 
combination  is  usually  called  the  NRL  Sum. 

Estimated  peaks  obtained  from  the  NRL  Sum  are  always 
positive  numbers,  but  are  to  be  interpreted  as  representing  peaks  that 
may  occur  in  either  direction. 


The  NRL  Sum  is  a  nonlinear  combination  and  must  always  be 
the  last  step  in  any  calculation.  In  particular,  if  the  NRL  Sum  is  used 
to  estimate  a  peak  acceleration  or  a  peak  force  from  a  combination  of 
modal  peaks,  that  estimated  peak  acceleration  or  force  cannot  be  used 
subsequently  to  find  stresses  or  deiormations. 

Estimating  peak  stresses  and  deformations 

Equations  18  and  19  were  especially  written  in  a  form  to  keep 
the  responses  of  the  structure  in  each  of  its  modes  separate  down  to 
the  point  of  calculating  stresses  and  deformations.  If  B  is  a  diagonal 
matrix  of  the  peak  accelerations  b  in  each  mode,  the  peak 
accelerations  of  the  structure  in  each  mode  are  obtained  by  scaling  the 
modeshape  matrix  according  to 

A  -  Q  B.  1 3  1 1 

The  peak  forces  on  the  structure  in  each  mode  are  given  by  the 
columns  of 

G  =  M  A  I32> 

The  stresses  and  deformations  for  each  mode  are  obtained  by  solving 
a  separate  static  problem  for  the  peak  loads  g  that  appear  in  each 
column  of  the  matrix  of  forces.  The  loads  mav  be  positive  or  negative, 
depending  on  the  signs  of  the  participation  factors  and  the 
modeshapes;  these  signs  must  be  preserved  through  Equations  31  and 
32. 


The  strain  energy  associated  with  a  peak  response  is 
proportional  to  the  square  of  the  strain  multiplied  by  the  elastic 
modulus  of  the  material,  or  proportional  to  the  square  of  the  stress 
divided  by  the  elastic  modulus.  Summing  the  squares  of  the  peak 
responses  thus  gives  a  number  proportional  to  the  total  energy  a 
particular  part  of  the  structure  must  accept  The  square  root  of  that 
sum  represents  the  single  deformation  or  stress  that  would  account 
for  the  same  energy 

The  NRL  Sum  provides  an  estimated  peak  response  that  is 
larger  than  the  energy-equivalent  peak,  but  not  larger  than  the  upper 
bound 


SELECTION  RULES  FOR  THE  NRL  SUM 


Peaks  to  be  included 


Peaks  from  the  following  sources  are  combined  into  a  single 
NRL  Sum: 

1  Stresses  and  deformations  from  normal  operation  of 
equipment. 

2  The  warping  deformations  and  stresses  from  differences  in 
motions  of  multiple  supports 

3  Stresses  and  deformations  from  some  but  not  usually  all  of 
the  normal  modes.  The  NRL  sura  overestimates  the  energies  for  the 
modes  that  are  included,  this  overestimate  can  be  considered  as  an 
allowance  for  additional  modes  of  the  structure  that  were  not  included 
in  the  sum 


1  For  simple  structures  that  can  be  modeled  with  six  or  fewer 
degrees  of  freedom,  the  three  modes  of  lowest  frequency  are  usually 
sufficient  for  determining  an  NRL  Sum. 

2.  Structures  of  moderate  complexity  (up  to  60  degrees  of 
freedom)  need  only  the  lower-frequency  modes  to  be  included  in  the 
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NRL  Sura.  About  half  as  many  modes  as  degrees  ot  Ireedora  is  an 
appropriate  combination  in  most  cases. 

3.  Complicated  models  with  very  many  degrees  of  freedom  las 
may  be  obtained  from  a  finite-element  model)  may  generate  a  modal 
thicket”  in  which  many  of  the  low-frequency  modes  represent 
vibrations  of  lightweight  parts  of  the  structure.  Here  it  may  be  better 
to  select  modes  in  accordance  with  modal  weight  rather  than 
frequency.  Modes  having  weights  at  least  2  percent  of  the  weight  of 
the  complete  model  are  most  likely  to  be  significant. 

Prestresses  and  bolted  icints 

Built-m  stresses,  including  prestresses  in  bolts,  are  not  included 
in  the  NRL  Sum  of  stresses  Bolted  joints  are  expected  to  respond  to 
tensile  'loads  by  decreasing  the  clamping  force  rather  than  by 
stretching  the  bolts  Special  checks  need  to  be  made,  however,  to  be 
sure  than  the  clamping  force  is  not  exceeded  by  a  peak  load  or  that 
shear  stress  from  a  transverse  load  does  not  combine  with  the  tensile 
stress  in  a  bolt  to  produce  an  excessive  value  of  maximum  normal 
stress. 

SUMMARY 

The  procedures  for  applying  the  theory  to  practical  calculations 
with  the  RIDDAM  can  be  summarized  in  eleven  steps: 

1.  Make  a  mathematical  model  of  the  structure. 

2.  Find  the  fixed-base  normal  modes  and  the  mass-normalized 
participation  factors  for  the  supports  of  the  model. 

3.  Add  the  participation  factors  for  each  mode  by  absolute 
value  and  square  the  sum  to  find  the  effective  mass  of  each  mode 

A  Refer  to  tables  and  formulas  showing  peak  accelerations  as  a 
function  of  the  weight,  frequency,  and  the  location  of  an  oscillator 
relative  to  the  point  of  impact  of  a  torpedo. 

3.  Combine  the  product  of  the  participation  factors  and  the 
tabulated  accelerations  for  each  support  to  find  the  peak  acceleration 
of  each  mode. 

6.  Scale  the  modeshapes  by  the  peak  accelerations  to  find  the 
peak  acceleration  of  each  point  on  the  structure  in  each  mode 


7  Multiply  the  peak  accelerations  of  the  structure  by  the  mass 
matrix  to  find  peak  forces  in  each  mode 

8.  Solve  the  static  problem  for  each  mode  to  find  peak  stresses 
and  deformations  at  critical  parts  of  the  structure. 

9.  Calculate  stresses  and  deformations  from  warping  of  the 
structure  caused  by  different  motions  of  redundant  supports.. 

10.  Use  the  NRL  Sum  to  combine  stresses  and  deformations 
from  operation,  warping,  and  selected  modes  at  critical  points  in  the 
structure. 

11.  Compare  the  estimated  peak  stresses  and  deformations 
with  allowable  values. 


SPECIAL  PROBLEMS 


DESIGN  SPECTRA 

The  accelerations  c(M,  P)  representing  the  peak  accelerations 
for  oscillators  having  masses  M  and  fixed-base  frequencies  P  can  be 
called  a  design  spectrum.  Each  point  in  a  design  spectrum  is  taken 
from  a  particular  time-history  response  of  a  massive  oscillator  That 
time  history  has  been  suppressed  for  simplicity. 

The  special  problems  described  here  arise  from  the  suppression 
of  the  time  histories  associated  with  the  design  spectra. 

REPEATED  FREQUENCIES 

Time -history  combinations 
Suppose 

ci  *  Pj2  zi  -  Pj2  u  1 55 1 

and 

C2  *  P22  Z2  -  P22  u  1 34) 

are  the  responses  of  two  oscillators  to  a  support  motion  urn  If  the 
two  frequencies  are  the  same.  Pj  =  P2  =  P  the  time  histones  will  be 
identical  and  the  peak  accelerations  will  occur  simultaneously 

As  the  accelerations  are  scaled  by  participation  factor  and  by 
raodeshape  and  then  used  to  calculate  forces,  stresses,  and 
deformations,  the  peaks  will  remain  coincident  m  time  but  may 
represent  responses  of  different  magnitude  occurring  in  either  the 
same  or  opposite  directions  The  combined  peak  is  given  by  the 
algebraic  sum  of  the  individual  peak  values 


1 
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NRL  Sum 

if  Modes  1  and  2  are  the  only  two  modes  of  the  structure,  the 
NRL  Sum  will  add  their  peak  stresses  and  deformations  by  absolute 
value  (absolute  value  of  the  larger  plus  square  root  of  the  square  of 
the  smaller).  This  is  an  appropriate  combination  if  both  responses  are 
in  the  same  direction,  but  it  may  be  a  very  large  overestimate  if  the 
modes  actually  produce  responses  in  the  opposite  directions 

If  the  structure  has  more  than  two  modes,  either  or  both  of 
Modes  1  and  2  may  be  part  of  the  square  root  of  the  sum  of  the 
squares  in  the  NRL  Sum.  Here  their  combined  response  will  be 
underestimated  if  they  respond  in  the  same  direction,  and  will  be 
overestimated  if  they  respond  in  opposite  directions. 

Special  rule 

Modes  with  identical  frequencies  should  have  their  peak 
responses  added  algebraically  rather  than  by  NRL  Sum. 

CLOSELY-SPACED  MODES 

DuhamdI's  tineara! 

Formal  solutions  for  Equations  33  and  34  can  be  written 


t 

cj  =  Pi2  J  v(t  I  cos[P i U  -  t  >)  dt 
0 


'33' 


and 


C2  =  P22  J  vi  t  1  co$!P2<t 
0 


1  i|  dt 


1  3b  1 


where 
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is  the  velocity  of  the  support. 


Time-historv  combinations 

When  the  two  accelerations  are  scaled  by  participation  t  actor 
modeshape,  mass,  and  stress  factor  they  will  produce  a  contribut: on 

t 

s(t)  -  Dj  J  v(t  i  cosfPi* t  -  t  ‘i  dt 
0 


t 

*  D2  f  vt  t  1  cos(P2‘t  -  t  • !  dt  r*5 ! 

0 

to  a  particular  stress,  where  Di  and  D2  are  combinations  0!  all  the 
relevant  factors  for  the  two  modes 

If  the  frequencies  are  Pi  -  P  -  E  and  P2  *  P  *  E.  with  an  average 
frequency  P  and  a  difference  of  2E.  the  stress  can  be  written 

1 

s't)  -  1 D 1  -  02'  |  vi  1  1  cosiE't  t  'I  cosIP't  t  '!  dt 
0 

t 

*  ' Di  -  D?'  |  v<t  1  sinlE't-t  '  siniP't-t  1)  dt 

0 


Approximation  for  closely-spaced  frequencies 

If  Et  is  less  than  tt/2,  the  factors  coslE'i  t  1]  and  smlF.'t  t  '!  can 
be  treated  as  positive -valued  windows  applied  to  the  velocity  vu  ' 
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sinlEf  1  -  cos(Et) 

with  average  values  — ^ —  and  - ^ -  over  the  range  of 

integration' from  0  to  t.  The  averages  can  be  considered  as  weighting 
factors  applied  to  the  velocity  of  the  support.  The  first  average  begins 
with  a  value  of  unity  at  the  beginning  of  the  response  and  decreases 
with  increasing  time,  while  the  second  begins  at  zero  and  increases. 


An  Et  of  IT/M  produces  an  average  of  0.900  for  the  first  window 
and  an  average  of  0  373  for  the  second.  For  frequency  differences  and 
times  up  to  this  limit  the  stress  mav  be  approximated  by 


si U  •  i D i  *  D21  J  vi t  t  coslP'l-t  'I  dt .  MO  1 

0 

to  reasonable  accuracy  That  is,  both  time  histones  approximate  the 
response  of  an  oscillator  at  the  average  frequency  P  and  their 
combined  stress  is  obtained  from  an  algebraic  sum  of  the  factors  Dj 
and  D2.  as  was  the  case  for  identical  frequencies. 

If  the  frequency  difference  2E  is  3  Hertz,  for  example.  Equation 
40  can  be  used  to  estimate  comf  med  peak  stresses  for  times  up  to  30 
milliseconds  alter  a  transient  motion  of  the  support  begins 


Convenience  is  served  u  the  responses  of  modes  having 
identical  or  closelv  spaced  irequencies  are  calculated  independently  in 
terms  of  their  weights  and  frequencies  Then  when  peak  responses 
are  to  be  combined  for  stresses  or  deformations  peaks  for  pairs  or 
clusters  o|  modes  having  neurlv  the  same  frequency  should  be 
combined  bv  algebraic  sum  to  produce  a  single  response.  That 
response  should  be  treated  as  the  response  of  a  single  mode  and 
combined  with  other  responses  bv  the  NRL  Sum 

Measurements  and  anaivses  of  the  responses  of  oscillators  to 
torpedo  impact  against  a  submarine  regularly  show  that  peak 
responses  occur  less  than  SO  milliseconds  after  the  initial  impact  The 
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peak  stresses  and  deformations  from  pairs  or  clusters  of  modes  that 
fall  within  a  5-Hertz  range  can  then  be  combined  algebraically  rather 
than  by  NRL  Sum. 

RULES  OF  THUMB 

Suppression  of  the  time  histories  is  seen  to  require  a  number  of 
fairly-arbitrary  rules  concerning  selection  of  modes  and  methods  of 
estimating  peak  values  from  combinations  of  responses. 

The  combination  rules  suggested  here  appear  to  be  based  on 
rational  criteria.  The  rules  were  chosen  mainly  to  avoid  obviously- 
mappropriate  results  and  are  subject  to  revision  and  reinterpretation 
in  light  of  statistical  results  from  experimental  data. 
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EXAMPLES 


A  TEXTBOOK  EXAMPLE 


Rigid  masses  weighing  325  pounds  and  200  pounds  are 
supported  by  springs  having  stiffnesses  of  10,000  and  30,000  pounds 
per  inch,  as  shown  in  Figure  1.  Each  mass  is  prevented  from  rotating 
by  a  set  of  roller  guides  that  are  not  considered  part  of  the  supporting 
structure.  The  masses  are  connected  by  a  thin  and  massless  beam  30 
inches  long.  The  cross-section  of  the  beam  has  an  area  moment  of 
inertia  of  1.25  in4  and  the  beam  is  made  of  a  material  with  elastic 
modulus  30  million  pounds  per  square  inch. 

The  structure  is  devised  to  illustrate  the  calculations;  it  is  not 
intended  to  represent  anything  realistic. 


For  coordinates  xpti  and  xj'ti  representing  displacements  of 
the  masses,  the  mass  matrix  is 


,r325  O: 

L  0  200j  P°unds- 


Moving  xj  bv  one  inch  and  holding  all  the  other  coordinates 
fixed  requires  a  force  of  10,000  pounds  to  stretch  the  spring  and  an 
additional  force  of 


(30E6 
<30  ji 


6.666  67  pounds 


to  bend  the  beam  as  a  double  cantilever  An  equal  and  opposite  force 
must  be  applied  to  X2  to  keep  the  other  end  of  the  beam  from  moving 
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Figure  1.  Textbook  structure.  Rigid  masses  are  labeled 
in  pounds.  Massless  springs  are  labeled  in  pounds  per 
inch.  The  masses  are  connected  by  a  massless  thin  beam 
30  inches  long  with  area  moment  1  25  inches  to  the 
fourth  power  and  elastic  modulus  30  million  pounds  per 
square  inch.  The  masses  are  on  fnctionless  guides  that 
prevent  them  from  rotating 
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Here  and  in  the  following  the  usual  computer  notation  of  E6  is  used, 
for  example,  to  indicate  multiplication  by  10  to  the  6th  power 

Moving  X2  while  the  other  coordinates  are  held  fixed  stretches 
its  spring  and  also  bends  the  beam.  The  forces  from  unit  displacement 
of  one  coordinate  at  a  time  can  be  assembled  into  the  stiffness  matrix 

["  26666.67  -16666.671 

K  “  L -16666.67  46666.67  J  Pounds  per  inch.  (431 


A  one-inch  displacement  of  ui  (with  everything  else  held  fixed) 
will  put  a  force  of  10,000  pounds  onto  xi-  A  similar  displacement  of 
U2  will  put  a  force  ol  30,000  pounds  onto  X2.  The  support-stiffness 
matrix  becomes 


The  three  matrices  M.  K.  and  F  allow  the  structure  to  be 
described  by  the  force-balance  equation  Ma  *  Kx  -  Fu.  shown  earlier 
as  Equation  1  in  the  section  on  Theory. 


The  next  step  in  the  analysis  is  to  find  a  matrix  Q  that  will  make 
QtMQ  be  a  unit  matrix  and  also  make  Q^KQ  be  a  diagonal  matrix,  as 
shown  in  Equations  3  and  4  This  can  be  done  most  easily  by  using  an 
available  computer  program,  such  as  in  the  example  shown  as  Case  1 
m  Appendix  B.  For  the  simple  structure  here,  the  modeshape  matrix  0 
can  be  calculated  by  hand,  using  the  formula  listed  in  Line  4  1 20  of  the 
computer  program  to  normalize  on  the  masses  and  the  formulas 
beginning  on  Line  3610  to  diagonalize  the  stillnesses.  The  modeshape 
matrix 
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ro.05 198516  -0.019351 14' 
1.0.02466795  0  06626833- 


i 


is  from  a  1 0 -digit  pocket  calculator.  A  direct  check  shows  that  it 

reduces 


-Tl__  ri.000  0.000] 

0  M  Q  «  L0.000  1.000J  =  1 


to  the  unit  matrix  and 


757.717  0.000' 

QTK  Q  =  L  0.000  257.663J  =  P 


to  a  diagonal  matrix. 

The  matrix  Q  can  be  used  to  represent  the  accelerations  and 
displacements  of  the  structure  in  terras  of  normal-mode  coordinates 
by  a  -  Qb  and  l  -  Qy.  as  shown  in  Equation  5 

Participation  factor 

The  participation-factor  matrix  is 


P-2  QTF 


r  9.007  1 2.82 1 1 


1-0.75 


7716. 


The  participation  factor  lets  the  equation  for  the  normal- mode 
coordinates  be  written  in  the  standard  form  b  *  p2y  =  pZRTu 
(Equation  7).  Each  element  of  Rtu  is  the  effective  displacement  of  the 
base  of  one  of  the  normal -mode  oscillators 

Modal  weights 


Add  along  the  rows  of  the  participation-! actor  matrix  hv 
absolute  value  and  square  the  result  to  obtain 

(9.007  *  1 2.82 1 '2  =  476  pounds  M9' 


v 

¥ 


for  Mode  1  and 


(0.751  +  7.716)2  =  72  pounds  150) 

for  Mode  2.  This  is  the  new  procedure  of  the  RIDDAM,  as  shown  by 
example  in  Equation  14. 

Mode  frequencies 

When  masses  are  in  pounds  and  stiffnesses  in  pounds  per  inch, 
frequencies  are  given  by  the  formulas 

y/386  (57.717)  =  23  8  Hertz  <51  ' 


and 


jjj'i 586  1257.668)  =  50.2  Hertz  <52< 

for  Modes  1  and  2.  The  normal-mode  equations  for  motionless 
supports,  b  *  P2y  >  0,  identify  the  diagonal  elements  of  Pz  as  squares 
of  the  fixed-base  natural  frequencies  of  the  modes.  The  factor  of  586 
is  the  acceleration  of  gravity  in  inches  per  second  per  second,  as 
needed  m  the  inch-pound-second  system  when  pounds  are  used  for 
both  masses  and  forces. 

Inputs 


The  weights  and  frequencies  of  the  modes  serve  as  entries  to 
tables  showing  the  peak  accelerations  of  the  modes.  As  an  example, 
suppose  that  the  tables  specify  a  peak  acceleration  of  26  times  the 
acceleration  of  gravity  <g'  for  a  weight  of  -176  pounds  and  a  frequence 
of  23.8  Hertz,  and  an  acceleration  of  1  19  g  for  72  pounds  at  50  2  Hertz 
These  are  values  at  the  point  of  impact  for  a  torpedo  At  a  point  on 
the  hull  30  inches  away  from  the  point  of  impact  the  response 
accelerations  would  be  1  1  70  g  and  53  55  g.  These  values  represent 
peak  accelerations  of  oscillators  of  mass  M  and  fixed-base  frequence  P 
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measured  during  tests  with  torpedo  impacts  or  caicuiated  from 
Equation  1 5  during  a  computer  simulation  of  an  impact. 

Peak  accelerations  of  the  modes 

If  the  torpedo  impacts  in  way  of  up  with  U2  located  30  inches 
away  from  the  impact,  the  peak  accelerations  of  the  modes  are 
obtained  by  multiplying  the  tabulated  responses  by  the  participation 
factors,  or 

9  007(26)+  12.821  (11.70)  =  384  1 9  g  1 3  3  ' 

for  Mode  1  and 

-0.751  (1  19) +  7.716  (53.55)  =  323.82  g  ( 54 ' 

for  Mode  2. 

Equations  53  and  54  correspond  to  applying  the  participation 
factors  to  the  predetermined  responses,  rather  than  applying  the 
participation  factors  to  the  support  motions  and  then  calculating  the 
responses.  The  two  procedures  give  equivalent  results  for  time- 
history  calculations,  as  explained  in  Equations  16  and  17  The  same 
combination  method  is  recommended  for  peak  values  in  Equation  30  if 
the  time  histories  have  been  suppressed. 

Peak  accelerations  of  the  masses 

Scale  the  shape  for  the  first  mode  (the  first  column  of  Qi  by 
384.19  and  the  shape  for  the  second  mode  by  323. 82  to  obtain 
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as  the  peak  accelerations  of  the  masses  in  each  mode  The  first 
column  of  A  shows  the  peak  accelerations  of  xi  and  x?  in  Mode  1  and 
the  second  column  shows  their  peak  responses  in  Mode  2.  Peak 
accelerations  are  shown  separately  for  each  mode,  as  recommended  in 
Equation  3 1 


Multiply  the  masses  (pounds)  and  the  accelerations  * g )  to  find 


M  .  r 6490  -20381 

M  A  "  '.1896  42902  pounds 


as  the  peak  forces  that  were  applied  to  each  mass  in  each  mode 


Each  mode  is  treated  as  a  static  problem  to  check  for  stresses 
and  deformations  at  critical  points  on  the  structure.  Suppose  that  the 
critical  points  are  the  bending  moments  at  each  end  of  the  beam  and 
the  tensile  loads  applied  to  bolts  restraining  the  springs  A 
redundant-structure  analysis  provides  deflections 


_  '0.3460  -0  0243" 

X  "  Lo  1642  0.0832J  incheS 
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for  each  mass  in  each  mode  when  the  forces  of  Equation  36  are 
applied  as  static  forces  to  the  structure.  The  deflections  produce 
forces  of  3460  pounds  at  ui  in  Mode  1  and  -243  pounds  in  Mode  2  At 
u?  the  forces  are  4926  pounds  in  Mode  1  and  2496  pounds  in  Mode  2 


Relative  motions  between  the  supports  will  distort  the  structure 
through  its  redundant  supports.  Suppose,  for  example,  that  the 
warping  formulas  specify  a  displacement  of  0  200  inch  at  the  impact 
point  and  a  displacement  of  0  090  inch  a  distance  of  30  inches  from 
the  impact  point  A  redundant-structure  anaivsis  shows  that  xi  will 
have  an  equilibrium  displacement  of  0  1431  inch  and  xz  a 
displacement  of  0.1090  inch  under  these  support  displacements  The 
resulting  tensions  in  the  hold-down  bolts  are  -369  pounds  at  Support 
1  and  an  equal  and  opposite  369  pounds  at  Support  2. 


Loads  on  the  supports  are  to  be  combined  by  NRL  Sum,  subject 
to  the  selection  rules  described  earlier  under  Application. '  The  forces 
from  warping  and  two  modes  (assuming  no  operating  stresses) 
combine  by  the  rule  of  the  largest  plus  square  root  of  the  sum  of  the 
squares  in  the  form 

3460  *  'J{- 243)2  -  (-369)2  =  4079  pounds  <58) 

as  an  estimate  of  the  peak  force  applied  to  the  bolts  at  uy  and 

4926  *  V(2496)2  ♦  (569)2  „  7486  pounds  (39) 

for  the  bolts  at  U2. 


The  peak  forces  of  Equations  58  and  59  are  interpreted  as 
applying  in  either  direction  to  the  hold-down  bolts.  The  bolts  need  to 
be  tight  enough  so  that  the  peak  force,  applied  in  tension,  will  not 
cause  the  joint  to  separate.  Similar  calculations  allow  peak  moments 
at  the  ends  of  the  beam  to  be  estimated  lor  comparison  with  allowable 
moments 

The  calculations  need  to  be  repeated  for  a  case  in  which  a 
torpedo  impacts  in  way  of  U2  rather  than  uy  The  most  severe 
conditions  can  be  taken  as  those  in  which  one  of  the  supports  for  an 
item  is  directly  inboard  of  the  point  of  impact  Then  a  check  of 
impacts  at  each  support  is  sufficient  to  evaluate  the  ability  of  the 
equipment  to  withstand  an  impact  at  any  point  near  it. 

AN  EXAMPLE  WITH  CROSS-AXIS  RESPONSE 

Structure 

A  point  mass  is  supported  by  two  equal  springs  battered  at 
angles  H  above  and  below  the  horizontal,  as  shown  in  Figure  2.  The 
mass  can  move  both  horizontally  t x i »  and  vertically  « X2 1  in  response  to 


i 


t-igure  i.  otruclure  with  cross-axis  response.  Mass  m  is 
supported  Dy  two  springs  of  stiffness  k  Each  spring 
slopes  at  an  angle  H  from  the  horizontal  The  dots 
represent  pinned  joints 
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the  horizontal  motions  of  its  two  supports.  This  is  another  textbook 
example  not  intended  to  represent  a  realistic  structure. 

In  this  and  in  the  following  two  examples  the  analysis  is 
presented  in  a  form  parallel  to  that  used  for  the  first  example,  but 
without  back-references  to  the  sections  of  the  report  on  theory  or 
application.  Please  refer  to  corresponding  sections  of  the  first  example 
for  appropriate  references. 

Matrices 


The  mass  matrix. 


M 


f  m  O' 

LO  mJ  • 


(60) 


represents  motion  of  the  mass  in  either  of  the  two  orthogonal 
directions.  The  stiffness  matrix 


f 2  k  cos2* Hi  0  ' 

K  '  L  0  2  k  sin2(H)  J 


(61 ) 


includes  factors  for  the  components  of  force  and  deflection  in  each 
spring  produced  by  small  displacements  of  the  mass  in  each  direction 
The  support-stiffness  matrix 

"  k  cos2*  H  i  k  cos2(  H )  ' 

^  “  L  -k  simH'i  cos(hTi  k  sin(H)  cosiHl  ' 


shows  the  horizontal  and  vertical  force  per  displacement  on  the  mass 
produced  by  small  horizontal  displacements  of  each  support. 

Normal  modes 


The  raodeshape  matrix 


v '.k rjrv  v  "J'  Jr f  rf  ^V^r'"/A'  "  'Av/"V-,«'-V-’.v.'>V/.  A  A  A  *”*  A  "<  J‘-  *’-  *"•  ■  ■.AAA<tAl''«.''«LA',‘'n.'>l 


normalizes  the  mass  matrix  to  a  unit  matrix,  gives  a  Irequencv- 
squared  matrix 

nTl,  A  i"2  (k/m)  cos2(H)  0  1 

9T K  9  *  L  0  2  ik/mi  sin2(H)J  'M' 


and  a  participation-factor  matrix 


p-2  qT  f 


7m/2 

( v  m/2 )  cot(H) 


_V  m/2 

i'/m/2i  cotiHi 


Modal  weights 

Mode  1  (horizontal  motion)  has  effective  mass  m,  as  obtained 
from  the  squared  sum  of  the  first  row  of  the  participation -factor 
matrix.  Mode  2  (vertical  motion)  has  effective  mass  m  cot2'H',  with 
the  square  of  the  cotangent  representing  the  effect  of  the  lever  arms 
that  act  to  magnify  or  diminish  the  effects  of  a  vertical  force  on  the 
horizontally-directed  supports  and  also  act  to  magnify  or  diminish  the 
motions  of  the  mass. 


Relative  horizontal  motions  between  the  two  supports  will 
simply  displace  the  mass  in  the  vertical  direction  without  producing 
any  stresses  or  deformations  (The  joints  at  the  supports  are  pinned  i 


The  simple  example  here  shows  that  the  matrix  formalism  is 
not  limited  to  unidirectional  responses  If  cross-axis  responses  are 
expected  from  unequal  motions  of  supports,  thev  should  be  included 
in  the  model  of  the  structure,  and  can  be  treated  on  the  same  basis  as 
on-axis  responses 


AN  EXAMPLE  WITH  ROTATIONAL  INERTIA 


Structure 

A  steel  plate  is  bridged  between  the  fianges  of  two  frames  to 
carry  a  weight  made  from  a  stack  of  steel  plates,  as  shown  in  Figure  5. 
Four  structures  Like  this  were  built,  installed  in  hulls,  and  tested  with 
torpedo  impacts,  The  one  analyzed  here  was  tested  in  1983  on  the 
SITV  (Submarine  Impact  Test  Vehicle)  There  it  served  as  simulated 
equipment  for  research  purposes. 


The  structure  was  especially  designed  to  be  symmetric  and 
nonredundant  with  two  simple  modes  of  vibration  In  one  mode  the 
weight  would  translate  as  a  rigid  body  and  bend  the  plate  as  a  beam 
The  stiffness  of  the  plate  <33  inches  iong,  22  inches  wide,  and  1  inch 
thick,  with  equal  loads  6.5  inches  from  each  end)  can  be  obtained  from 
handbook  formulas  as 


12  (28E6)  22 
(6.5)2  (3(33)  - 


1)3/12 
4(6. 5)1 


99,724  pounds  per  inch 


The  stack  of  plates  (859  pounds  '  and  the  center  20  inches  of  the 
flexible  plate  (124  pounds)  both  participate  in  the  bouncing  motion, 
while  the  6.5  inches  at  each  end  of  the  flexible  plate  connect  from  the 
fixed  support  to  the  moving  load  'sav  one-third  of  81  pounds'  A 
rough  estimate  o.  the  total  mass  in  the  bouncing  mode  is  thus  1010 
pounds.  The  frequency  should  be 


1  ,  386  M99724 

2n  ;  loio 


4 A  Hertz 


Analysis  of  measurements  from  the  SITV  test  showed  that  the 
bouncing  mode  had  a  fixed-base  frequency  of  44  Hertz 
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WEIGHT 


Forward 


Starboard 


Figure  3  Structure  with  rotational  inert. na  The  weight 
(a  stack  of  steel  plates)  is  Doited  to  a  flexible  plate  that 
bridges  between  the  flanges  of  two  successive  ring 
stiffeners  of  a  submarine  hull  Half-round  spacers  (not 
illustrated)  provided  nearly  moment-free  connections  at 
the  ends  of  the  flexible  plate  and  at  the  connections 
between  the  plate  and  the  weight 


Rocking 

If  the  weight  rotates  about  an  axis  through  the  center  ol'  the 
plate  and  parallel  to  its  ends,  each  half  of  the  plate  will  bend  as  a 
simply-supported  beam  with  a  point  load.  Handbook  formulas  give 
the  stiffness  of  each  half  as 


3  (28E6)  (22  (1)3/12]  16.5 
(6.5  )2  (16.5  -  6.5)2 


-  601,420  pounds  per  inch.  (68 1 


Two  such  loads  are  each  10  inches  from  the  center  of  the  plate,  and  a 
small  rotation  produces  deflections  of  10  inches  per  radian.  The 
angular  stiffness  is  thus  120.284  million  pound-inches  per  radian. 

The  mass  moment  of  inertia  for  rotations  about  an  axis  through 
the  center  of  the  plate  can  be  obtained  by  starting  with  the  moment 
about  the  center  of  gravity  of  the  weight  (a  block  25  5  by  7  by  17 
inches), 


1 25.5)2  *  (7)2 

859 - — -  =  50,055  pound-in2  1 6 9  t 

Transfer  to  a  parallel  axis  through  the  middle  plane  of  the  plate  i9 
inches  from  the  center  of  gravity  of  the  weight)  bv 

859  (9 12  =  2 1 .475  pound-m2,  (70 1 

and  add  an  estimate  for  the  center  part  of  the  plate 

„  1 20  *  ( ]  i2  ^ 

124  ^ -  -  4, 1 42  pound-in2.  (71  i 

to  obtain  a  total  moment  of  inertia  of  75,682  pound-in2  Frequencv  is 
then 


J_  ,  ,386  ( 1 20.284E6 ) 
2 tt  v  75682 


125  Hertz. 


(72) 


Data  from  the  test  showed  a  fixed-base  frequency  of  137  Hertz  for  the 
rocking  mode 


Matrices 


The  mass  matrix. 


r  1  o  i  o 


M  “  i  0  75682J 

contains  a  mass  in  pounds  and  a  mass  moment  of  inertia  m  pound 
inches  squared.  The  stiffness  matrix, 


,r  199724 
L  0 


20.284E6J- 


has  elements  in  pounds  per  inch  and  in  pound -inches  per  radian.  The 
support-stiffness  matrix 

f  99862  99862  ] 

F  *  L3.645E6  -3.645E6J  '  1 

shows  that  a  one-inch  displacement  of  both  supports  would  apply  a 
force  of  199,724  pounds  to  the  translational  coordinate. 
Displacements  of  +1  and  -1  inch  would  cause  a  rotation  of  1/16  5 
radians  and  apply  a  moment  of  1 20  2S4E6/16.5  pound-inches  to  the 
rotational  coordinate. 


The  modeshape  matrix 

;3146584E-3  0  ^ 

Q  ’  L  0  3  63499E-3-  76 

normalizes  the  mass  matrix  to  a  unit  matrix  and  gives  a  frequency 
squared  matrix 


VVv«  *  *.  *.  * » ^  * .  • .  * .  i . .. .  . .  a  . . . . -  -  - 


QtKQ 


f  1 97.746 

L  0.000 


0  000  l 
1589. 334j  = 


that  gives  the  frequencies  of  AA  and  125  Hertz. 
Participation  factor 


The  participation-factor  matrix  is 


P-2  QT  F 


'15.8903  15  8903'  T 

L  8.3365  -8.3365 J  =  R 


177) 


Add  along  the  rows  by  absolute  value  and  square  the  sums  to  find 
that  Mode  1  (bouncing)  has  effective  weight  1010  pounds  and  Mode  2 
(rocking)  has  effective  weight  278  pounds. 

Responses 

From  this  point  onward,  the  analysis  proceeds  in  the  usual 
fashion.  The  1DDAM  inputs  at  each  support  for  the  first  mode  are 
obtained  from  its  weight  (1010  pounds),  its  frequency  <AA  Hertz),  and 
the  location  of  each  support  relative  to  the  point  of  impact  of  the 
torpedo.  Multiply  each  input  by  its  factor  from  the  first  row  of  the 
participation-factor  matrix  and  sum  to  find  the  acceleration  of  the 
mode.  Scale  the  first  column  of  the  raodeshape  matrix  by  that 
acceleration  to  obtain  the  acceleration  in  the  bouncing  mode  of 
response  and  continue  to  calculate  forces,  stresses,  and  deformations 

Calculations  for  the  second  mode  will  be  similar  except  that  its 
accelerations  will  be  rotational  accelerations  in  units  of  386  radians 
per  second  squared.  Multiplv  by  the  mass  moment  of  inertia  in 
pound-inches  squared  to  obtain  moment  in  pound-inches,  and 
continue  to  calculate  stresses  and  deformations  from  the  applied 
moment. 

Finally,  the  NRL  Sum  is  used  to  combine  stresses  from  the  two 
modes  with  warping  stresses  and  operating  stresses  to  estimate  peak 
stresses  at  critical  points.  The  particular  structure  of  the  example  had 
neither  warping  nor  operating  stresses  Checks  were  made  for  tensile 
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loads  in  the  bolts,  stresses  in  the  flexible  plate,  and  stresses  in  the 
welded  brackets. 

Interpretation 

The  example  here  shows  that  the  formalism  applies  as  well  to 
rotational  responses  as  it  does  to  translations,  with  the  rotations  given 
in  radians,  stiffnesses  represented  by  moments  per  radian,  and 
masses  represented  by  mass  moments  of  inertia. 

Alternate  models  using  translations 

Figure  4  shows  four  ways  in  which  the  two  coordinates  (one 
translation  and  one  rotation)  for  the  example  could  be  replaced  by  a 
pair  of  translational  coordinates.  In  each  case  the  total  mass  of  1010 
pounds  is  represented  by  point  masses  located  so  as  to  reproduce  the 
mass  moment  of  inertia  about  an  axis  through  the  center  of  the  plate 
The  case  marked  la'  is  of  special  interest,  since  it  shows  that 
placing  tw’o  weights  of  130  pounds  each  directly  over  the  support 
points  not  only  gives  the  correct  mass  moment  of  inertia  but  provides 
a  total  of  278  pounds  to  match  the  effective  mass  of  the  rocking  mode 
Each  of  the  configurations  in  Figure  4  can  be  analyzed  to 
produce  results  identical  to  the  results  obtained  using  a  rotational 
coordinate. 

A  REAL  EXAMPLE 

Source 


Figure  5  is  adapted  from  a  sketch  in  Appendix  E  of  a  report  of 
engineering  analyses  of  submarine  equipment  done  by  a  Navy 
shipyard.1  The  figure  shows  a  4 -cubic-foot  airflask  supported 
between  two  hull  stiffeners  bv  four  beams  built  of  quarter-inch  steel 


*M  A  HattamerandM  W  Yargus  Analysis  of  MSW  Cooling  Pump  &.  Piping  and 
Miscellaneous  Air  Flasks  using  Dynamic  Analysis  Method  for  Impact 
Bremerton.  Washington  Puget  Sound  Naval  Shipyard  Engineering 
Document  449  (22  April  1 077 ) 
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Figure  4,  Representing  rotation 
translotions.  Rectangles  are  do 
masses  in  pounds  Heavu  sirale 
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Figure  5.  Sket.cn  or’  moael  tor  an  a:r  flask  The  flask  is 
mounted  between  frames  Du  four  beams  made  of 


quarter-inch  plate 
shown 


Cross-sect’ons  of  the  beams  are 
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plate.  Calculations  from  the  analysis  are  summarized  here  with  some 
changes  in  notation  and  revisions  in  procedures. 

Coordinates 

The  airflask  was  taken  to  have  three  degrees  of  freedom, 
represented  by  an  athwartship  translation  of  the  center  of  gravity 
< x  1 ).  rotation  about  a  vertical  axis  (xj).  and  rotation  about  a  fore-and- 
aft  ax::  1x3). 

Mass  matrix 


Masses  were  represented  in  units  of  386  pounds  (slmches),  as 
required  by  the  computer  program  that  was  to  be  used  to  find  normal 
modes.  The  mass  for  the  translational  motion  was  1.17  slmches.  The 
mass  moment  of  inertia  about  the  vertical  axis  (along  the  axis  of  the 
flask)  was  taken  as  that  for  a  thin  circular  shell  with  radius  9  inches 
or 


1.17  (9)2  -  95  lb-in-s2.  (79) 

The  shell  was  taken  as  42  inches  long  for  its  moment  of  inertia  about 
the  fore-and-aft  axis,  or 

6(9)2  *  4f42)2 

1.17 - — -  -  117  1 2 1  I2  =  2  !  9  lb -in  ~s2.  1 80  ' 

The  mass  matrix  was  thus 
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Beams 

The  two  beams  connecting  to  Frame  90  were  angles  made  of 
quarter-inch  plate  4  inches  long,  3  inches  wicie.  and  6  inches  deep 
The  cross-sectional  area  was  2.28  square  inches  and  the  area  moment 


about  the  neutral  axis  was  9  4  in4  Their  bending  stiffness  as 
cantilevers  fixed  to  the  frame  and  pinned  to  the  flask  would  be 


-30E6(9.4)  ....  . 

o  — — ^ —  =  lp.2  million  Ib/in. 

Shear  stiffness  (neglecting  any  shear  factor)  would  be 


2.25(1 1E6) 


6.2  million  lb/in. 


giving  a  combined  stiffness  of 


132(6.2) 
,3  2  *  6.2 


=  4.2  million  ib/in 


for  each  beam. 

The  two  beams  connected  to  Frame  89  were  12  inches  long.  3 
inches  wide,  and  11  inches  deep,  with  area  3  5  square  inches  and  an 
area  moment  of  46.4  in4  The  combined  stiffness  in  bending  and  shear 


2.413.2) 
2  4  *  3  2 


4  million  lb/in. 


All  four  beams  act  against  an  athwartship  displacement  of  the 
flask,  producing  a  stiffness  of  112  million  pounds  per  inch  for 
Coordinate  1  The  beams  are  attached  7  inches  forward  and  7  inches 
aft  of  the  centerline  of  the  flask,  so  that  when  the  flask  is  rotated 
about  a  vertical  axis  they  apply  a  moment  of  548.8  million  pound- 
inches  per  radian  to  Coordinate  2  Moment  arms  for  rotation  about  a 
fore-and-aft  axis  are  12  inches  above  and  below  the  center  of  the 
flask  for  a  stiffness  of  1612  8  million  pound-inches  per  radian  in 
Coordinate  3 


if"  V 


In  addition,  when  the  flask  is  displaced  to  starboard  (out  of  the 
paper  in  Figure  5)  the  beams  will  apply  a  moment  about  the  vertical 
axis  of 


(-4.2-4.2M7)  +  (-1.4-1.4H-7)  =  -39.2  million  lb-in/in  (86) 

that  must  be  reacted  by  an  equal  and  opposite  moment  to  keep  the 
flask  from  rotating.  Similarly,  rotation  about  the  vertical  axis  will 
produce  a  net  force  per  angle  of  39.2  million  pounds  per  radian  to  port 
that  must  be  reacted  by  an  equal  force  to  starboard  to  keep  the  flask 
from  translating.  The  stiffness  matrix  needs  to  be  written  as 


T11.2E6 
K  =  I39.2E6 

L  o 


39.2E6 
348. 8E6 
0 


1 6 1 2.8E6J 


A  displacement  of  Frame  90  in  the  starboard  direction  at  the 
upper  support  (uj)  will  apply  4.2  million  pounds  of  force  per  inch  of 
displacement  to  Coordinate  1  The  force  acts  through  a  moment  arm 
of  7  inches  to  apply  29  4  pound-inches  of  moment  per  inch  of 
displacement  to  Coordinate  2.  and  acts  through  a  moment  arm  of  12 
inches  to  apply  50.4  million  pound-inches  per  inch  to  Coordinate  3. 

Similar  calculations  for  the  other  three  supports  lead  to  a 
support-stiffness  matrix 

j  4.2E6  4.2E6  1  4E6  1.4E6 

F  =  !  29.4E6  29.4E6  -9.8E6  -9  8E6  :  (88) 

l50.4E6  -50.4E6  16.8E6  -16.8E6, 


Normal  modes 

Data  were  entered  into  Computer  Program  RGGGl  using  the 
special  format  it  required.  The  program  returned  three  modes  with 
natural  frequencies  298,  424,  and  536  Hertz  It  provided  a 
modeshape  matrix. 


sVW’.'* 
V  "V  -  v  v 


rw*3r*jr*y 


r  i.  o.  1.  I 

QS  =  !  0.185  0.  -0.066  \  (89) 

[  0.  1.  0.  j 

that  had  been  scaled  to  make  the  largest  element  in  each  column  have 
a  value  of  unity. 

Program  RGGG1  did  not  display  the  scale  factors  that  it  used. 
They  were  recalculated  by  hand  in  a  tabulation  scheme  that  was 
equivalent  to 

r4.42  0  0  ' 

S2  =  (QS)t  M  QS  »  0  219  0  .  i90) 

L  0  o  1  -58 j 

The  original  analysis  continued  with  rounded  calculations  to  remove 
the  squares  of  the  scale  factors  from  individual  results  as  thev 
appeared.  The  following  calculations  depart  from  that  method  to 
remove  the  scale  factors  from  the  modeshape  matrix  directly  by- 
dividing  each  column  by  the  square  root  of  its  element  of  S2.  to  obtain 

ro.476  0  000  0.796  ' 

Q  *  i  0.088  0.000  -0  053  '  (91  ' 

L0  000  0  068  0.000  j 


Alternate  calculation 


A  recalculation  of  the  normal  modes  was  oiade  using  a  1 0 -digit 
pocket  calculator  to  obtain  a  mass-normalized  modeshape  matrix 
directly  as 


,r  0  483 
Q  -  -0  087 
'  0.000 


0.000  -  0  788 
0.000  0  054 

0  068  0.000 


1 92  1 


Differences  in  the  third  significant  figure  between  this  modeshape 
matrix  and  the  rescaled  matrix  from  RGGGi  appear  to  result  from 
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rounding  of  the  scaled  modeshapes  in  Equation  89..  Frequencies  of  the 
modes  are  given  as  298.  432.  and  848  Hertz  from  the  revised  analysis. 

The  signs  are  reversed  for  rotation  about  the  vertical  axis  in 
Equation  92  because  the  data  were  entered  into  RGGG1  using  a  left- 
handed  coordinate  system.  (Direction  x  pointed  aft.  y  upward,  and  z  to 
starboard.) 

The  modeshape  matrix  from  Equation  92  has  been  used  in  the 
following  calculations. 

Participation  factor 

The  participation  factor  is 

r -0.1 56  -0.186  0.438  0.438  1 

rT  .  p-2  qT  p  «;  0.462  -0.462  0.154  -0.154;.  1 9  3 ) 

:  0.412  0.412  0.049  0.049  j 

Modal  weights 

Sum  each  row  of  the  participation-factor  matrix  by  absolute 
value  and  square  the  sums  to  find  modal  masses  in  the  units  of  386 
pounds  used  in  the  mass  matrix.  Multiply  bv  386  to  find  modal 
weights  of  550  pounds  for  Mode  1.  587  pounds  for  Mode  2.  and  328 
pounds  for  Mode  3. 

IDDAM  inputs 

One  of  the  IDDAM  tables,  used  for  the  sake  of  an  example, 
shows  that,  at  the  point  of  impact,  a  weight  of  550  pounds  mounted  at 
298  Hertz  would  respond  with  a  peak  acceleration  of  1264  g.  A 
weight  of  587  pounds  at  432  Hertz  would  respond  to  1431  g.  and  328 
pounds  at  548  Hertz  would  respond  to  2297  g. 

One  frame  away  from  the  impact  point  the  responses  would  be 
0  45  times  as  large,  and  at  a  distance  of  24  inches  around  a  frame 
from  the  impact  point  the  responses  would  be  0.60  times  as  large 


If  a  torpedo  were  to  strike  in  way  of  the  upper  support  at  the 
after  frame  (support  u j  l.  the  IDDAM  inputs  for  the  four  supports, 
multiplied  by  the  participation  factors,  predict  peak  accelerations  for 


the  three  modes  of 

["84  0 

0  ; 

B  -  1  0  304 

o  ; 

(94) 
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The  scaled  modeshape 

matrix 

is 

r  -40 

j 

0 

1259: 

A  *  Q  B  •  i  -7 

0 
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(98  1 

L  0 

21 

0  j 

with  accelerations  in  g  for  the  translational  coordinate  (first  row  of  the 
matrix)  and  angular  accelerations  in  multiples  of  386  radians  per 
second  per  second  for  the  rotational  coordinates  (second  and  third 
rows). 

Forces  and  moments 
Multiply 

'18  0  868  ' 

G  =  0.386  M  A  -  -269  0  3136  i96i 

L  0  1739  0  J 


to  obtain  peak  forces  (thousands  of  pounds)  applied  to  Coordinate  1 
and  peak  moments  (thousands  of  pound-inches'  applied  to 
Coordinates  2  and  3  in  each  of  the  three  modes  (The  factor  of  0  386  is 
needed  to  return  the  mass  matrix  to  convenient  units  i 


The  force  and  moment  trorn  Mode  1  place  a  load  of  -5  thousand 
pounds  on  each  of  the  two  after  supports,  and  a  load  of  14  thousand 
pounds  on  each  of  the  forward  supports. 

The  moment  from  Mode  2  applies  a  load  of  36  thousand  pounds 
to  each  of  the  upper  supports  and  -36  thousand  pounds  to  each  of  the 
lower  supports. 

Mode  3  puts  loads  of  254  thousand  pounds  on  each  after 
support  and  30  thousand  pounds  on  each  forward  one. 

Warning 

Suppose  that  the  warping  formula  gives  a  displacement  0.200 
inch  at  the  impact  point  ( u i ),  0.120  inch  at  U2,  0.090  inch  at  U3,  and 
0.054  inch  at  u^.  The  static  formula  for  the  motion  of  the  airflask  can 
be  solved  most  easily  by  writing  it  in  the  form 

x  -  K  1  F  u  =  iQ  P  2  qT)  F  u  =  Q  Rt  u  (97) 

that  avoids  inverting  the  stiffness  matrix.  The  equilibrium  positions 
of  the  coordinates  are 

'0.1  16000; 

x  -  0.006286  :  (98  * 

L0  002875  - 

The  displacement  (translation  and  two  rotations)  at  support  u; 
is  0  194502  inch,  deflecting  the  beam  there  by  0  008498  inch  and 
producing  a  force  of  23  thousand  pounds  Similar  calculations  show 
balancing  forces  of  -23  thousand  pounds  at  U2.  -23  thousand  pounds 
at  U3,  and  *23  thousand  pounds  at  U4. 

Stress  analysis 


The  outermost  fibers  of  the  beams  at  Frame  90  were  4  inches 
from  the  neutral  axis  and  would  have  stress 


(99) 


<  9^  =  !  702 

psi  in  bending  per  pound  of  force  applied  at  the  end  of  the  beam. 
Shear  stress  would  be 


25-  -  0.400 


(100) 


psi  per  pound. 

The  beams  at  Frame  89  had  their  outermost  fibers  6  inches 
from  the  neutral  axis,  with  stress  1733  psj  per  pound  there  in 
bending,  and  with  shear  stress  0.286  psi  per  pound. 


The  upper  aft  beam  f u  1 )  has  loads,  in  thousands  of  pounds,  of  -3 
•.Mode  1),  36  (Mode  2).  234  (Mode  3).  and  23  'warping'  Bending 
stresses  are  -9,  62.  433.  and  39  thousand  pounds  per  square  inch  The 
NRL  Sum  of  the  bending  stresses  is 


433  +  V(-9)2  *  ( 6 2 )2  ♦  (59)2  -  507  thousand  psi 


MOb 


Shear  stresses  are  -2.  14,  102.  and  9  thousand  pounds  per 
square  inch,  with  an  NRL  Sura  of  1 19  thousand  psi 


Check 

The  calculation  might  as  well  stop  here  'as  it  did  in  the  original' 
because  there  is  no  chance  that  the  support  at  the  impacted  frame 
could  withstand  the  bending  stress  of  507  thousand  pounds  per 
square  inch  or  the  shear  stress  of  1  19  thousand  pounds  per  square 
inch  The  bending  stress  is  more  than  ten  times  the  yield  stress  for  a 
medium  steel  and  the  shear  stress  is  more  than  four  times  yield  in 
shear.  A  redesign  of  the  supports  lor  the  airflask  is  obviously  needed 


APPENDIX  A 

A  SHORT  INTRODUCTION  TO  MATRIX  ALGEBRA 


MATRICES 

A  computer  statement  similar  to 

DIMENSION  M(4,3>  (A-H 

will  reserve  20  memory  spaces  for  numbers  that  can  thereafter  be 
referenced  as  M(J,L).  with  J  an  integer  in  the  range  from  1  to  4  and  L 
an  integer  in  the  range  from  1  to  5  The  numbers  are  usually  thought 
of  as  arranged  in  a  rectangular  array  with  four  rows  and  five  columns, 
designated  simply  by  the  letter  M: 

"Ml  1.1)  Mil, 2)  MC  1.3)  MI  1.4)  M(  1.5')“ 

MI2.1)  MI 2.2)  MI2.3)  MI2.4)  MI2.5) 

M  ’  MI 3.1 )  MI 3 .2 )  MI3.3)  MI. 3.4)  MI3.5)  iA'^ 

LM<4.  1 )  MI 4.2)  MI 4.3 1  MI 4.4)  MI4.5H 

SOME  DEFINITIONS 

The  matrix  with  dimension  MI 4.3  *  is  called  a  rectangular  maim 
because  it  has  a  different  number  of  rows  than  columns.  A  matrix 
with  dimension  MI  4.4)  would  be  called  a  square  matrix  The 
elements  M(J,J)  are  called  the  diagonal  elements ,  they  would  be  the 
elements  M>  1.1),  MI 2,2),  MI 3.3).  and  Ml 4.4)  for  the  matrix  here 

The  sum  of  the  diagonal  elements  is  called  the  trace  of  the 
matrix;  the  trace  of  the  matrix  here  would  be  MI  1,1 )  *  Ml 2.2'  -  VI '3.3' 
*  MI 4.4 »  A  symmetric  matrix  is  one  m  which  every  MiJ.Ll  is  equal  to 
the  corresponding  MIL.J)  -  that  is.  Mil. 2'  =  M ( 2 . 1  •  Mil. 3'  =  M« 3.1  i. 
and  so  on  for  all  the  pairs 

SOME  OPERATIONS 

If  matrix  Q  has  elements  QlJ.L).  the  matrix  with  elements  QiL.J) 
is  called  the  transpose  of  Q  designated  by  Transposition 

corresponds  to  interchanging  row?  and  columns  in  the  arrav.  ii  the 


transposed  matrix  is  equal  to  the  original  matrix  -  Q),  then  the 
matrix  must  have  been  both  square  and  symmetric. 

Matrices  of  the  same  dimension  are  added  or  subtracted  by 
adding  or  subtracting  the  corresponding  elements.  Matrices  of 
different  dimensions  cannot  be  added  to  or  subtracted  from  one 
another. 

MATRIX  MULTIPLICATION 

The  product  of  two  matrices  is  written  as  if  it  were  ordinary 
multiplication.  That  is, 

Z  =  M  Q  (A-3) 

is  the  product  of  M  and  Q.  The  operation  is  not  an  element-by- 
element  multiplication,  however,  but  involves  multiplying  each 
element  in  a  row  of  the  first  matrix  by  the  corresponding  element  in  a 
column  of  the  second  matrix  and  summing  the  products  Each  element 
of  Z  is  generated  as 

Z(J.N)  -  M(J,1 )  Q(  1  ,N)  -  M(J.2)Q(2.N)  ♦  M(j.3)Q<3.N' 

♦  M(J,4)  QH.N)  *  ....  (A-4) 

The  equation  can  be  abbreviated  as 

Z(J.N)  -  MiJ.Ll  QlL.N  i.  '  A  -  5  ' 

where  it  is  understood  (the  summation  convention  i  that  the  product 
is  to  be  summed  over  all  possible  values  of  the  repeated  index  L. 

The  product  matrix  Z  will  have  as  many  rows  as  there  were 
rows  m  M.  and  as  many  columns  as  there  were  columns  in  Q 
Moreover.  M  musi  have  exactlv  as  many  columns  as  Q  has  rows,  so 
that  the  index  L  can  be  assigned  to  both  and  run  through  a  single 
range  of  values 

Because  of  the  rule  for  multiplication,  it  is  not  surprising  that 
the  matrix  product  QM  is  a  different  matrix  from  the  product  MQ 


MORE  DEFINITIONS 


A  diagonal  matrix  is  one  in  which  ail  of  the  elements  other  than 
the  diagonal  elements  have  values  of  zero.  A  unit  matrix  (usually 
designated  by  I)  is  a  diagonal  matrix  in  which  each  diagonal  element 
has  a  value  of  unity. 

The  inverse  of  a  matrix  K  is  the  matrix  K*1  that  makes 

K'1  K  =  I  (A -6) 

be  a  unit  matrix.  Not  all  matrices  have  inverses.  The  definition  of  an 
inverse  is  reciprocal;  K  is  also  the  inverse  of  K1  so  that 

UK1  =  I  (A -7) 

If  the  transpose  of  a  matrix  is  also  its  inverse,  the  matrix  is  said 
to  be  unitary  That  is. 

UTU  -  I  i a i 

defines  U  as  a  unitary  matrix.  A  simple  example  is 
fcos(H)  -smlH.H 

U  -  .sin(H)  cos' H ) J-  <A'<51 

which  is  unitary  for  any  angle  H. 

VECTORS 

A  matrix  with  only  one  column  is  called  a  column  vector  If  a 
matrix  X  has  dimension  X'5,1 ).  for  example,  it  is  a  column  vector  with 
five  elements.  It  is  usual  to  suppress  the  unit  designation  and  show 
its  dimension  as  X( 5 J  Column  vectors  are  designated  by  lower-case 
boldface  letters,  such  as  x,  in  the  present  report 

A  matrix  with  only  one  row  is  called  a  row  vector  It  is 
dimensioned  X»  1.5 for  example,  and  is  designated  here  using  the 
symbol  for  transpose,  x^  Again  the  unit  designation  is  usuallv 


suppressed  to  show  a  row  vector  as  dimension  X(5>.  although  this  has 
the  defect  of  not  distinguishing  it  from  a  column  vector. 

All  of  the  rules  for  adding,  subtracting,  and  multiplying 
matrices  apply  to  vectors  provided  the  unit  dimension  is  understood. 

SOME  HELPFUL  RELATIONS 

The  following  rules  can  be  verified  by  working  out  some  simple 
examples.  First,  the  transpose  of  a  matrix  product  is  obtained  from 
the  transposes  of  the  factors  taken  in  the  reverse  order,  or 

(MQ)T  -  QTMT  i A- 10) 

Similarly,  the  inverse  of  a  product  involves  inverting  and 
interchanging  the  factors,  as 

(MQ)'1  -  Q'1  M-l.  •  A- 1 1  i 

If  a  matrix  is  symmetric,  us  inverse  is  also  symmetric 
However,  the  product  of  two  symmetric  matrices,  such  as  M'1  and  K. 
is  not  a  symmetric  matrix,  since 

(M*1  K)T  -  Kt  M-T  -  K  M  1  ( A  - 1 2 1 

is  different  from  the  original  product  M  1  K 

Both  the  mass  and  stiffness  matrices  M  and  K  for  a  structure 
are  necessarily  symmetric  If  a  stiffness  matrix  were  not  symmetric, 
for  example,  the  structure  could  be  run  as  a  perpetual-motion 
machine  by  moving  it  around  a  cycle  in  which  the  smaller  off-diagonal 
element  was  worked  during  the  compression  phase  and  the  larger 
offdiagonal  element  was  used  to  provide  power  during  the  return 
stroke 

Matrix  elements  for  a  structure  can  often  bv  generated  or 
checked  bv  using  the  useful  relations  that  x"  K  x/2  is  the  potential  u>r 
strain)  energy  and  vT  M  v/2  is  the  kinetic  energy,  where  x  is  a 
column  vector  of  displacements  and  v  a  column  vector  of  velocities 


USEFULNESS  OF  MATRIX  ANALYSIS 

The  rules  for  addition,  subtraction,  and  multiplication  of 
matrices  were  chosen  specifically  so  that  matrix  equations  very 
similar  to  scalar  equations,  such  as 


&  i  =  F  u. 


(A-I3) 


wouid  expand  into  precisely  the  set  of  coupled  differential  equations 
needed  to  represent  a  structure  with  multiple  coordinates.  The  great 
compaction  of  the  notation  and  the  resultant  improvement  in 
comprehension  are  exactly  what  is  needed  in  analyzing  complicated 
structures. 
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APPENDIX  B 

A  COMPUTER  PROGRAM  FOR  RIDDAM 


PROGRAM 


1000  REn  THIS  PROGRfln  IS  URITTEN  IN  CORE  BPS  1C 


1100  REH  I  DENT  I FV  THE  PROGRfln 
1110  PRINT  "877F7  RIOOAfi" 

1120  PRINT  "  ACCEPTS  R  0 1  AGONAL  HASS  HATRIX  IN" 

1130  PRINT  "  POUNDS  OR  POUND  INCHES  SQURRED." 

1110  PRINT  “  STIFFNESSES  IN  POUNDS  PER  INCH  OR" 

1150  PRINT  M  POUND- INCHES  PER  RADIAN." 

1160  PRINT  "R  L  BORT,  NAUAL  RESEARCH  LAB,  UASHINGTON,  DC" 
1170  PRINT  " JULV  15,  1937." 

1180  PRINT 


1200  REn  D I  HENS  I  ON  THE  UARIRBLES 


:,M  M"c,io),  kg: 

1220  Din  0(10,10),  P(1Q, 10),  R( 1 0, 1 0) 

1230  Din  FI (10) ,  U(10),  0(10),  F2( 1 0) , B( 1 0) 

1210  Din  AGO, 10), G(  10, 10) 

1300  REn  CASE  NUnBER  OR  TEST  CASE 

1310  PRINT  "ENTER  A  NUHBER  FOR  VOUR  CASE,  OR" 

1320  PRINT  "ENTER  0  (ZERO)  TO  RUN  A  TEST  CASE." 

1330  INPUT  T 

1310  IF  T >0  THEN  1110 

1350  PRINT  "CASE  0  IS  A  1000-POUND  BAR  SUPPORTED" 

1360  PRINT  “  INDEPENDENTLY  AT  EACH  END.  IT  IS  HOOELED  AS" 
1370  PRINT  “N1  EQUAL  HASSES  UlTH  Ml  *  1  EQUAL  SPRINGS." 

1330  PRINT 
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MOO  REH  ASSIGN  DEGREES  OF  FREEDOM 

1410  PRINT  "ENTER  DEGREES  OF  FREEOGO  (1  TO  10) 

1 120  INPUT  N1 

1430  IF  HI < 1  THEN  1410 

1440  IF  HI >10  THEN  1410 

1500  REfl  CLERR  SOME  flflTR I CES 

1510  FOR  J*1  TO  N1 

1520  FOR  L=1  TO  N1 

1530  LET  n(J,L)30 

1540  LET  K( J, L)=0 

1550  LET  Q(J,L)=G 

1560  NEXT  L 

1570  NEXT  J 

1500  IF  T>0  THEN  1810 

1600  REM  TEST  CASE,  NRSS  RNO  STIFFNESS 
1610  LET  N2=2 
1620  FOR  J-l  TO  N1 

^  1  “  ^  V*  '  -  •  '  '  l  •  ^  « 

.  O«j0  Ul  II  V  S< ,  O  J  ~  i  JUU/  v  il  I  I  i 

1640  LET  K( J, J)=2*20000*(N1  +1 ) 

1650  NEXT  J 

1660  FOR  J=1  TO  Ml -1 

1670  LET  K(J,J+l)=-20000*(N1+1) 

1680  LET  K(  J+1 , J )=K( J , J+1 ) 

1690  NEXT  J 

1700  REH  TEST  CASE,  SUPPORT  STIFFNESS 

1710  FOR  J*1  TO  N1 

1720  LET  F(J, 1 )=0 

1730  LET  F(J,2)=0 

1740  NEXT  J 

1750  LET  F(1,1)»20000*(N1+1) 

1760  LET  F ( N 1 , 2)=20000*(N1  +1 ) 

1770  GOTO  2410 
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1800  REH  ASSIGN  NUMBER  OF  SUPPORTS 

1810  PRINT  "ENTER  NUMBER  OF  SUPPORTS  (1  TO  10) 

1820  INPUT  N2 

1830  IF  N2< 1  THEN  1810 

1840  IF  N2>10  THEN  1810 

2000  REH  ENTER  THE  MATRIX  ELEMENTS 


$ 

;v 

a* 


1 


2100  REM  DIAGONAL  MASS  MATRIX  M 
2110  FOR  J-1  TO  N1 
2120  PRINT  "ENTER  MASS"; J 
2130  INPUT  M( J, J) 

2110  IF  FI ( J , J ) >0  THEN  2170 

2150  PRINT  “MASSES  MUST  BE  PQSITIUE.  TRY  AGAIN." 
2160  GOTO  2120 
2170  NEXT  J 

2200  REM  SYMMETRIC  STIFFNESS  MATRIX  K 
2210  FOR  J-1  TO  N1 
2220  FOR  L-J  TO  N1 

2230  PRINT  "STIFFNESS  MATRIX  ROU",.J;  "COLUMN"  ;L 
2240  INPUT  K(J,L) 

2250  LET  K(L,J)«K(J,L) 

2260  NEXT  L 
2270  NEXT  J 

2300  REM  SUPPORT-STIFFNESS  MATRIX  F 
2310  FOR  J-1  TO  N1 
2320  FOR  L-1  TO  N2 

2330  PRINT  “STIFFNESS  COORDINATE" ; -J;  "TO  SUPPORT" ; L 
2340  INPUT  F(J,L) 

2350  NEXT  L 
2360  NEXT  J 
2370  PRINT 
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2400  REN  LABEL  THE  CASE 
2410  PRINT  “THIS  IS  CASE";T 
2420  PRINT 

2500  REN  ECHO  THE  HASS  NATRIX 
2510  PRINT  "HASS  HATR I X  n  =“ 

2520  FOR  J-1  TO  HI 
2530  FOR  L-1  TO  HI 
2540  PRINT  n ( J , L ) , 

2550  NEXT  L 
2560  PRINT 
2570  NEXT  J 
2580  PRINT 

2600  REH  ECHO  THE  STIFFNESS  HATR IX 
2610  PRINT  "STIFFNESS  nATRIX  K  «" 

2620  FOR  J-1  TO  N1 
2630  FOR  L-1  TO  N1 
2640  PRINT  K(J,L), 

2650  NEXT  L 
2660  PRINT 
2670  NEXT  J 
2680  PRINT 

2700  REH  ECHO  THE  SUPPORT-STIFFNESS  NATRIX 
2710  PRINT  "SUPPORT-STIFFNESS  NRTRIX  F  «“ 
2720  FOR  J-1  TO  N1 
2730  FOR  L-1  TO  N2 
2740  PRINT  F(J,L), 

2750  NEXT  L 
2760  PRINT 
2770  NEXT  J 

2790  PRINT 
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3000  REM  CHANGE  OR  CORRECT  HfiTR I X  ELEHEHTS 
3100  REH  CHECK  THE  MATRICES 

3110  PRINT  "PLEASE  CHECK  THE  flfiTRICES.  IF  THEV  RRE  GOOD" 
3120  PRINT  "ENTER  0.  IF  V0U  UANT  TO  CHRNGE  RNV,  ENTER  1 . " 
3130  INPUT  J 
3140  IF  J=0  THEN  4110 

3200  REH  CHANGE  THE  NASS  MATRIX  n 

3210  PRINT  "NASS  TO  BE  CHANGED  OR  0  IF  NO  CHANGE" 

3220  INPUT  J 

3230  IF  J<1  THEN  3310 

3240  PRINT  H( J, J); " IS  OLD  nA$S";J;\  ENTER  NEU  UALUE ; " 
3250  INPUT  n(J,J) 

3260  IF  n( J, J) >0  THEN  3210 

3270  PRINT  “MASSES  MUST  BE  POSIT IUE  NUMBERS !“ 

3230  GOTO  3240 

3300  REM  CHANGE  THE  STIFFNESS  MATRIX  K 

3310  PRINT  "ENTER  RGU,  COLUMN  OF  ELEMENT  TO  BE" 

3320  PRINT  "CHANGED  IN  K.  ENTER  0,0  IF  NO  CHRNGE." 

3330  INPUT  J,L 
3340  IF  J=Q  THEN  3410 

3350  PRINT  K(J,L); “ IS  0L0";J;L;".  ENTER  NEU  UALUE:" 

3360  INPUT  K( J, L) 

3370  LET  K ( L , J ) =K ( J , L ) 

3380  GOTO  3310 
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3400  REH  CHARGE  THE  SUPPORT-ST IFFMESS  MRTRIX  F 
3410  PRINT  "ENTER  ROU,  COLUMN  OF  ELEMENT  TO  BE" 

3420  PRINT  "CHANGED  IN  F.  ENTER  0,0  IF  NO  CHRNGE. " 

3430  INPUT  J,L 

3440  IF  J*0  THEN  2410 

3450  PRINT  F(J,L); “ IS  0L0“;J;L;“.  ENTER  NEU  URLUE: " 
3460  INPUT  F(J,L) 

3470  GOTO  3410 


4000  REN  START  NORMAL-MOOE  RNfiLVS  I S 

4100  REH  NORMALIZE  ON  THE  MASS  MATRIX,  QT*M*Q  = 

4110  FOR  >1  TO  N1 

4120  LET  Q(J,JH/SQR(M(J,J)) 

4130  NEXT  J 


4200  REM  SVMMETRIC  OVNAMIC  MATRIX,  OT*K*Q 
4210  FOR  J=1  TO  N1 
4220  FOR  L*1  TO  N1 

4230  LET  P ( J ,  L ) aQ ( J , J ) *K ( J , L ) *Q ( L , L ) 

4240  NEXT  L 
4250  NEXT  J 


4300  REM  SHOU  THE  OVNAMIC  MATRIX 

4310  PRINT  “SVMMETRIC  OVNAMIC  MATRIX  QT*K*Q  =" 


4320  FOR  >1  TO  N1 
4330  FOR  La1  TO  N1 
4340  PRINT  P( J, L) , 
4350  NEXT  L 
4360  PRINT 
4370  NEXT  J 
4380  PRINT 


5000  REfl  DIAGONAL  1 2E  THE  OYNAMIC  HATRIX 


5100  REH  SET  COUNTER  AND  RUN  THROUGH  OFFDIRGONAL  ELEMENTS 
5110  PRINT  “CALCULATING. . 

5120  LET  E=0 
5130  FOR  J*1  TO  N1-1 
5110  FOR  L=J+1  TO  N1 

5200  REM  SORT  COLUMNS  BY  INCREASING  FREQUENCY 
5210  IF  P(J, J)<=P(L,L)  THEN  5110 
5220  FOR  J1*1  TO  N1 
5230  LET  2=Q(J1 ,  J) 

5210  LET  Q(J1,J)-Q(J1,L) 

5250  LET  Q(J1,L)-Z 
5260  LET  Z«P(J1,J) 

5270  LET  P(J1,J)-P(J1,L) 

5280  LET  P(J1,L)«2 
5290  NEXT  J1 

5300  REH  AND  INTERCHANGE  ROUS  OF  P  RS  UELL 
5310  FOR  J1-1  TO  HI 
5320  LET  Z=P(J,J1) 

5330  LET  P(J,J1)«P(l,J1) 

5310  LET  P(L, J1  )*2 
5350  NEXT  J1 

5100  REM  SKIP  SMALL  URLUES 

5110  IF  RBS(P(J,L)+P(L, J))< .0!  THEN  5920 

5500  REM  DON'T  DIUIOE  BY  ZERO 
5510  IF  P(J,J)<P(L,L)  THEN  5610 
5520  LET  C=SOR( .5) 

5530  LET  S-C 
5510  GOTO  5710 


:,T.'^Vily>>vv',v.vxviv.v.\».-.«.,A’  V  •_«■ 


5600  REH  JACOB  I  ROTATION 

5610  LET  H-.5*ftTH((P(JlL)*P(LiJ))/(P(J,J)-P(LJL))) 
5620  LET  C»COS(H) 

5630  LET  S-SIH(H) 

5700  REN  UPOATE  THE  COLUNNS  OF  THE  NATRICES 

5710  FOR  J1-1  TO  HI 

5720  LET  Z»C*Q(J1,J)+S*Q(J1,L) 

5730  LET  Q(J1,L)=C*Q(J1,L)-S*Q(J1,J) 

5740  LET  g(J1,J)»2 

5750  LET  2=C*P( J1 , J)+S*P( J1 , L) 

5760  LET  P ( J 1 ,L)=C*P(J1 , L)-S*P(J1 , J) 

5770  LET  P( J1 , J)*Z 
5780  NEXT  J1 


5800  REN  AND  ALSO  THE  ROUS  OF  P 

5810  FOR  J1-1  TO  HI 

5820  LET  ZaC*P( J, J1 )+S*P(L, J1 ) 

5830  LET  P(L,J1)=C*P(L,J1)-S*P(J,J1) 

5840  LET  P(J,J1)«Z 
5850  NEXT  J1 

5900  REN  STEP  THE  COUNTER  AND  END  THE  LOOP 
5910  LET  £*■£«■  1 
5920  NEXT  L 
5930  NEXT  J 

6000  REN  SHOU  AND  TELL  THE  NORDAL  NODES 
6100  REN  IS  P  DIAGONAL  VET? 

6110  PRINT  "ZAPPED" ;E; "OF  THE  OFFD I  AGONAL  ELEHENTS." 
6120  IF  E>0  THEN  5120 
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6200  REH  SHOU  flOOESHflPES 

6210  PRINT  "flOOESHRPE  flATR I X  0  =" 

6220  FOR  J-1  TO  N1 
6230  FOR  L»1  TO  HI 
6240  PRINT  Q(J,L), 

6250  NEXT  L 
6260  PRINT 
6270  NEXT  J 
6280  PRINT 

6300  REM  SHOU  FREQUENCY-SQUARED  fIRTR  I X 

631 Q  PRINT  “FREQUENCY -SQUARED  HftTR I X  QT*K*Q 

6320  FOR  J-1  TO  N1 

6330  FOR  1*1  TO  N1 

6340  PRINT  P(J,L), 

6350  NEXT  L 
6360  PRINT 
6370  NEXT  J 
6380  PRINT 

6400  REH  CHECK  FOR  LOOSE  PARTS 

6410  FOR  J-1  TO  N1 

6420  IF  P( J, J)< . 0  THEN  6440 

6430  LET  P(J, J)*-.01 

6440  NEXT  J 

6500  REH  PARTICIPATION-FACTOR  T1ATRIX  R 

6510  FOR  J-1  TO  N1 

6520  FOR  L-1  TO  N2 

6530  LET  R( J,L)=G 

6540  FOR  J1*1  TO  HI 

6550  LET  R ( J ,  L ) *R ( J ,  L ) +Q ( J 1 ,  J )  *F ( J 1 ,  L ) 

6560  NEXT  J1 

6570  LET  R(J,L)=R(J,L)/P(J, J) 

6580  NEXT  L 
6590  NEXT  J 
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|  6600  REfl  SHOU  THE  PftRT I C I PfiT I  OH  FfiCTOR 

j  6610  PRINT  "PRRTICtPRTlOH  FACTOR  R  =“ 

I  6620  FOR  J=1  TO  N1 

6630  FOR  L=1  TO  N2 
6610  PRINT  R(J,L)( 
f  6650  NEXT  L 

j  6660  PRINT 

!  6670  NEXT  J 

6680  PRINT 


6700  REH  CHECK  FOR  NE6RTIUE  EIGEHUflLUES 

6710  FOR  J-1  TO  N1 

6720  IF  P(J,J)>0  THEN  6770 

6730  IF  P(J,J)>-.011  THEN  6760 

6710  PRINT  "STRUCTURE  COLLAPSES  IN  nODE";J 

6750  PRINT  "THERE  IS  R  HI  STAKE  IN  VOUR  STIFFNESS  MATRIX, 

6760  LET  P( J, J)=0 

6770  NEXT  J 


^  W  V-V  WW  ?X?X  *1'  ■J^V* '.’, ', 


K 


7000  REH  FREQUENCIES,  UEIGHTS,  RND  INPUTS 

7100  REN  REAL  OR  FAKE  INPUTS 
7110  IF  T=G  THEN  7410 

7120  PRINT  "ENTER  0  TO  GET  FAKE  INPUTS  (EXRtlPLE  GNLV) , 
7130  PRINT  “OR  ENTER  1  TO  SUPPLV  INPUTS  FROH  TABLES." 

7 HO  INPUT  J 

7150  IF  J-0  THEN  7410 

7200  REH  ENTER  INPUTS 
7210  FOR  J«1  TO  N1 

7220  PRINT  "INPUT  FOR" ; U(J) ; "POUNDS  AND" ; FI (J) ; "HERTZ ; 
7230  INPUT  D(J) 

7240  NEXT  J 

7300  REfl  ENTER  FAROE  FACTORS 
7310  FOR  J-1  TO  N2 

7320  PRINT  "FRAHE  FACTOR  FOR  SUPPORT" ;J 
7330  INPUT  F2( J) 

7340  NEXT  J 
7350  GOTO  7610 

7400  REH  INPUTS  (FAKE)  FOR  TEST  CASE 
7410  FOR  J»1  TO  N1 

7420  LET  0( J)«. 3*P(J,J)/$Qfi( 1+. 000001 *U(j)*P(j, j) ) 

7430  NEXT  J 

7500  REfl  FRAHE  FACTORS  (FAKE)  FOR  TEST  CASE 

7510  LET  F20H 

7520  FOR  J=2  TO  N2 

7530  LET  F2(J)«.5*F2(J-1) 

7540  NEXT  J 
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7600  REH  SHOU  THE  HODES 

7610  PRINT  ,,nOOE,,/'FRE0UENCV“, '‘UE I GHT"/' INPUT- 
7620  PRINT  , "(HERTZ)"; “(POUNDS)", “(G)" 

7630  FOR  J-1  TO  HI 

7640  PRINT  J,  FI (J) ,U(J),D(J) 

7650  NEXT  J 
7660  PRINT 

7700  REH  SHOU  THE  FRfiflE  FACTORS 

7710  PRINT  "FRANE  FACTORS  (RATIO  TO  II1PACT  POINT) 

7720  FOR  J-1  TO  N2 

7730  PRINT  F2(J), 

7740  NEXT  J 
7750  PRINT 
7760  PRINT 

8000  REH  PEAK  RESPONSES 

8100  REH  CALCULATE  ACCELERATIONS  OF  NODES 

8110  FOR  J-1  TO  HI 

8120  LET  B ( J ) =0 

8130  FOR  1=1  TO  N2 

8 MO  LET  B( J)=B( J)+R( J, L)*F2(L)*0( J) 

8150  NEXT  L 
8160  NEXT  J 

8200  REN  CONUERT  TO  ACCELERATIONS  OF  HASSES 

8210  FOR  J-1  TO  N1 

8220  FOR  L=1  TO  Nt 

8230  LET  A( J, L)=Q(J,L)*B(L) 

8240  NEXT  L 
8250  NEXT  J 
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8300  REN  FI  NO  PERK  FORCES  IN  ERCH  NODE 

8310  FOR  J=1  TO  N1 

8320  FOR  L=1  TO  NT 

8330  LET  G(J,L)=N(J,J)*fi(J,L) 

8310  NEXT  L 
8350  NEXT  J 

8100  REN  SHOU  THE  RCCELERRTIONS 

8110  PRINT  “PERK  RCCELERRTIONS  (G)  IN  ERCH  NODE" 

8120  FOR  J-1  TO  N1 

8130  FOR  L-1  TO  N1 

8110  PRINT  R(J,L), 

8150  NEXT  L 
8160  PRINT 
8170  NEXT  J 
8180  PRINT 

8500  REN  SHOU  THE  FORCES 

8510  PRINT  "PERK  FORCES  (POUNDS)  IN  ERCH  NODE- 

8520  FOR  J-1  TO  N1 

8530  FOR  L-1  TO  N1 

8510  PRINT  G(J,L), 

3550  NEXT  L 
8560  PRINT 
8570  NEXT  J 
8580  PRINT 

9000  REN  UINDUP 


9100  REft  UHRT  3HRLL.  UE  00  HE^TT 

9110  PRINT  “THIS  IS  ALL  THIS  PROGRRN  DOES.  CHOOSE: 


9120  PRINT  -1 
9130  PRINT  "2 
9110  PRINT  "3 
9150  PRINT  "1 
9160  PRINT  "5 
9170  PRINT  ”6 
9180  PRINT  "9 


NEU  INPRCT  POINT" 

NEU  INPUT  TABLES" 

CHANGE  THE  SUPPORTS" 

CHANGE  THE  STIFFNESS  NRTRiX" 
CHANGE  THE  NRSS  NRTRIX" 

NEU  STRUCTURE- 
END  THIS  RUN" 


9200  REO  CHOICE 
9210  INPUT  J 
9220  IF  J=1  THEN  7310 

9230  IF  J=2  THEN  7210 

9240  IF  J=3  THEN  3410 

9250  IF  J=4  THEN  3310 

9260  IF  J=5  THEN  3210 

9270  IF  J=6  THEN  1310 

9280  IF  J=9  THEN  9910 

9290  GOTO  9110 

9900  REH  END 

9910  PRINT  “THANK  VOU  FOR  RUNNING  THIS  PROGRfitl . 
9920  PRINT  “FINISHED- 
9930  END 

COMMENTS 


Documentation 

The  program  is  self-documented  by  remarks.  It  uses  the 
same  procedures  and  most  of  the  same  symbol?  as  are 
described  in  the  text  of  this  report. 


Core  basic 


The  program  is  written  in  primitive  BASIC  language 
without  any  of  the  enhancements  that  are  available  m  various 
dialects  of  BASIC.  As  such,  it  is  long  and  slow,  but  should  run  as 
is  on  any  computer  that  has  a  BASIC  compiler  or  interpreter 
The  program  is  also  wasteful  of  memory  since,  for  clarity  it 
does  not  use  any  shortcuts  and  defines  new  variables  lor  each 
step  in  the  calculation. 
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SAMPLE  OUTPUT 

Test  case.  Case  0 

The  program  is  loaded  with  an  example  consisting  of  a 
uniform  bar  with  weight  1000  pounds  and  axial  stiffness 
20,000  pounds  per  inch.  The  bar  can  be  represented  by  from 
one  to  ten  lumped  masses  evenly  spaced  along  a  massless 
spring  A  simple  analytical  solution  is  available  for  the  fixed- 
base  frequencies  and  modeshapes  of  the  string  of  equal 
masses. - 

Following  is  part  of  the  output  from  Case  0  The  bar  was 
modeled  using  four  stations,  each  representing  the  motion  of  a 
200-pound  mass  lump.  The  remaining  mass  was  lumped  with 
the  supports  at  each  end  of  the  bar 

The  test  case  uses  a  simple  formula  to  generate  1DDAM- 
itke  inputs.  The  inputs  are  for  sake  of  the  example  only  and  are 
not  intended  to  apply  to  any  particular  situation. 

THIS  IS  CASE  0 


HASS  hRTRIX  n 

= 

200 

0 

0 

0 

0 

200 

0 

0 

0 

0 

200 

0 

0 

0 

0 

200 

STIFFNESS  HfiTR 

IX  K  • 

200000 

-100000 

0 

0 

-100000 

200000 

-100000 

0 

0 

-100000 

200000 

-’00000 

0 

0 

-100000 

200000 

-Fan  Y  Chen  On  Modeling  and  Direct  Solution  of  Certain  Free  Vibration 
Systems.'  Journal  of  Sound  and  Vibration  Volume  M  Number  I  U'-Vli 
pages  "57  to  79 
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SUPPORT-STIFFNESS  HRTR I X  F  = 

1 00000  0 

0  0 

0  0 

0  100000 

svnnETRic  DVNflmc  nflTR  i  x  qt*k*q  = 
1000  -500  0 

-500  1000  -500 

0  -500  1000 

0  0  -500 

CALCULATING. . . . 

ZAPPED  6  OF  THE  OFFD I  AGONAL  ELEflENTS 
ZAPPED  6  OF  THE  OFFD I  AGONAL  ELEMENTS 
ZAPPED  5  OF  THE  OFFD I  AGONAL  ELEMENTS 
ZAPPED  0  OF  THE  OFFD I  AGONAL  ELEMENTS 


0 

0 

•500 

1000 


£ 

MODESHAPE  MATRIX  Q  = 

2.623614E-G2 

-4.25328E-02 

4 . 253267E-02 

-2. 628634E-02 

s 

4.253229E-02 

-2 . 628693E-02 

-2 . 623676E-02 

4 . 25324E-02 

» 

4.25328E-02 

2.628614E-02 

-2.628635E-02 

-4.253267E-02 

t 

2 . 628698E-02 

4.253229E-02 

4 . 253241 E-02 

2 . 628677E-02 

fc 

$ 

FREQUENCY-SQUARED  MATRIX  QT*K*Q  = 

J 

i 

190.933 

4.936953E-03 

2.26589E-05 

5 . 789629E-06 

i 

4.912365E-03 

690.983 

-5 . 086626E-06 

-2 . 625585E-05 

3 . 227984E-05 

4.204654E-06 

1309.017 

-2 . 479627E-03 

L* 

C 

c 

-7 . 09665E-06 

2.624094E-05 

-2 . 425924E-03 

1809.017 

V 

>. 

1 

"H 

PARTICIPATION 

FACTOR  R  * 

13.7635 

13.76404 

M 

M 

-6.155405 

6 . 1 55331 

■; 

3.249208 

3.249187 

!■ 

r 

-1.453074 

1  .453097 

HO 


4 


3 


£ 

jf 

«n 

f. 


a 


4 


;* 

:? 

■Si 


•X 

V 

I 


% 

V, 

V 


nooE 

FREQUENCV 

UEIGHT 

I NPUT 

(HERTZ) 

(POUNDS) 

(G) 

i 

13,21263 

757.771 1 

53.55082 

2 

32.19539 

151 .5512 

197.2251 

3 

113.1323 

12.22911 

382.2812 

1 

132.9919 

8.115327 

538.6061 

FRflHE  FACTORS 

1 

(RATIO  TO 
.5 

IHPRCT  POINT) 

PERK  RCCELERRT 

IONS  (G)  1 

IN  EACH  NODE 

29.0617 

25.81773 

79.21529 

10.28613 

17.02323 

15.95639 

-13,97652 

-16.61339 

17.02331 

-15.95588 

-18,97571 

16.61319 

29.06262 

-25.81712 

79.2118 

-10.2863 

PERK  FORCES  (POUNDS)  IN 

EACH  NODE 

5812.31 

5163.516 

15819.06 

2057.226 

9101.655 

3191 .279 

-9795.301 

-3328.678 

9101.768 

-3191.177 

-9795.117 

3328.693 

5812.521 

-5163.131 

15813.96 

-2057.259 

Following  is  output  from  the  problem  described  as 
Textbook  Example  in  the  text. 


THIS  IS  CASE  I 


rmss  nPTR  i  x  n  = 

325  0 

0  200 


STIFFNESS  HATRIX  K  » 
26666.67  -16666.67 

-16666.67  *6666.67 


ami  wjt  "jwv 


"Wf.w^v 


SUPPORT-ST I FFMESS  MATRIX  F  = 

10000  0 

0  30000 

SYMMETRIC  DYNAMIC  nfiTR  I X  QT*K*Q  = 
82.0513  -65.37206 

-65.37206  233.3331 


CRLCULAT I NG . . . . 

ZAPPED  1  OF  THE  OFFD I  AGONAL  ELEflENTS. 
ZAPPED  0  OF  THE  0FF0 1  AGONAL  ELEflENTS . 


H00ESHAPE  MATRIX  Q  = 

5. 198516E-02  -1 .9351 1 1E-02 

2.166796E-02  6.626833E-02 


FREQUENCY-SQUARED-  MATRIX  QT*K*Q  = 
57.71698  0 

1 . 907319E-06  257.6677 


PARTICIPATION  FACTOR  R  = 
9.006909  12.82185 
-.7510112  7.715557 


MODE  FREQUENCY  UEIGHT 

(HERTZ)  (POUNDS) 

1  23.75558  176.1916 

2  50.19307  71.68277 

FRAME  FACTORS  (RATIO  TO  IMPACT  POINT) 
1  .15 


INPUT 

(G) 

26 

119 


PEAK  ACCELERATIONS  (G)  IN  EACH  MOOE 
19.97215  -6.265853 

9.177313  21.15753 


s: 


*  .  •  .  • 


PEAK  P0ACE6  (POUNDS)  IN  £ACM  NODE 

6491.047  -2036.402 

1895.463  4291.506 


